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Abstract 



Q ' We review recent developments concerning the weak decay of A hypernuclei. 

2 . New studies covering the mesonic decay channel as well as recent models for 

the non-mesonic one are discussed and compared with experimental data. The 
puzzle of the neutron- to proton- induced decay ratio, Tn/Tp, is addressed in 
^ ■ connection to the two-nucleon induced decay channel and proposals for more 

efficient experimental analyses of this ratio are made. 
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1 Introduction 



Hypernuclear physics is reaching the stage of a mature science. Many theoretical and 
experimental efforts have been devoted to the subject and a number of comprehensive 
reviews are already available |I|,0,^,^,|^,|^,|^,|^,|^,[T^,[TI|]. The present review covers 
only one aspect of the investigations in the field, that of the weak decay of A hyper- 
nuclei. This topic has proved to be a rich one in the interface between particle and 
nuclear physics: both information on particle properties unaccesible with ordinary 
elementary reactions, as well as nuclear properties of hypernuclei, complementary to 
those from nuclear structure, have been obtained. Many of the interesting results 
from hypernuclear decay have been extensively reported in the reviews quoted above, 
hence the present one will be devoted essentially to the new achievements available 
from 1990 on. 

The early experiments to produce hypernuclei using emulsions or the (i^~,7r^) 
reaction at CERN and BNL have given room to many new ones in a large number of 
experimental facilities: S and A hypernuclei have been produced at KEK using the 
{K~ , it) reaction with stopped kaons. The same reaction but with in-flight kaons has 
been studied at BNL. More recently, the {'n'^,K~^) reaction at BNL and KEK has 
proved to be very efficient in populating deeply bound A states and, when studied 
away from the forward scattering angles, has allowed the measurement of asymmetries 
in the particles emitted from the decay of polarized hypernuclei. Another source of 
information has been obtained from LEAR, where hypernuclei and their decay are 
studied using p absorption followed by delayed fission. A similar technique but with 
protons on nuclei has been successfully developed at COSY which has allowed the 
measurement of lifetimes of heavy hypernuclei. The TJNAF laboratory has also joined 
these efforts and the (e, e'K^) reaction is now ready to produce hypernuclei with a 
much better energy resolution. Double strange AA and S hypernuclei have also been 
obtained from the {K~, K^) at BNL, opening new possibilities to the understanding 
of the A A and HiV interactions as well as to unravel the existence of the H particle. 
Finally, the FINUDA project at DA$NE will soon use the tagged and slow K~ 
from the $ decay into K^K'^ to provide data at a rate significantly higher than was 
available in the past. 

The theoretical developments have run in parallel to the experiments. The mesonic 
decay of A hypernuclei (A — >■ vriV) has confirmed the strong sensitivity to the pion 
nucleus optical potential discovered in the past. The total mesonic decay rate is sig- 
nificantly enhanced due to the pion interaction in the nucleus which comes, basically, 
from the attractive p-wave part of the optical potential. However, exclusive reactions 
to a final closed shell nucleus select basically the repulsive s-wave part and lead to a 
reduced partial decay rate. The simultaneous study of inclusive and exclusive decay 
reactions can thus be a good source of information on the pion nucleus interaction, 
complementary to the one obtained from pionic atoms and low energy pion nucleus 
scattering. The mesonic decay of light hypernuclei has provided evidence of the strong 
repulsion of the AA^ interaction at short distances, a property that follows naturally 
from available quark models of the strong YN interaction. 

The non-mesonic decay mode (AA^ — > NN) has received much attention and 
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detailed models which go beyond the original one pion exchange mechanism are now 
available, either including the exchange of several mesons or working explicitly in 
terms of the quark degrees of freedom. One of the motivations for developing these 
models was the large discrepancies between theory and experiment in the ratio F^/Fp 
of neutron- to proton-stimulated non-mesonic decay, where the one pion model seems 
to provide an order of magnitude smaller ratio than the data. The question is not yet 
settled but one of the things that has been recently found is that the consideration 
of the two-nucleon-induced decay mode {ANN NNN) enlarges the error bars 
from present experimental analyses. The two-nucleon-induced non-mesonic decay of 
A hypernuclei was considered with the hope that it could solve the F„/Fp puzzle. 
While it has been seen that this is not the case, the new decay remains, however, as 
a channel which must be taken into account in any attempt of a precise analysis of 
the experimental data. Theoretical studies of the neutron and proton spectra from 
the decay of A hypernuclei are now available and open new doors for a more reliable 
determination of the F„/Fp ratio than was possible in the past. 

In this review we concentrate only on the weak decays of A hypernuclei. The 
progress in E hypernuclei has been scarce both experimentally and theoretically, 
however, the narrow S states claimed by old spectra have not been observed in a 
recent experiment carried out at BNL with better statistics. Another interesting 
experiment has been the clean measurement of the |He hypernucleus at BNL ||12|| , 
giving stronger grounds for the existence and properties of this state than the previous 
findings at KEK using stopped kaons. Also, the availability of AA and S hypernuclei 
has stimulated many theoretical studies aimed at obtaining information about the 
AA and SN interaction. A comprehensive recent review on these latter issues can be 



found in ref. 11 
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2 Weak mesonic decay of hyper nuclei 



2.1 Derivation of the A width in nuclei 

We make here a formal derivation of the width of A states bound in nuclei. In the 
first place we evaluate the A width for a A particle moving through infinite nuclear 
matter. The width for finite nuclei is obtained from there using the local density 
approximation (LDA). In a second step we follow a direct approach to the evaluation 
of the A mesonic width in finite nuclei and compare the methods. Anticipating results 
we find that the nuclear matter plus LDA is a good tool to evaluate A decay widths 
in nuclei. For the non-mesonic channel, where the energy carried by the emitted 
nucleons is large, the LDA predicts decay rates within 5% of the finite nuclei results. 
The mesonic width is more sensitive to nuclear shell effects given the little phase space 
available for the reaction. In medium hypernuclei around C, assuming one uses the 
correct Q value of the reaction, one may expect the LDA results to lie within 20-30 % 
of the correct answer. In lighter hypernuclei the mesonic width approaches gradually 
the free one and the uncertainties naturally diminish. In heavier hypernuclei, where 
the mesonic rate is reduced in several orders of magnitude, the errors are larger. In 
this case the LDA can only account qualitatively for the mesonic width, however it 
provides the essential features of the reaction, such as the reduction in the rate with 
increasing mass of the hypernucleus, as well as the dramatic increase of the rate due 
to the pion interaction with the nucleus. 

We also discuss here two ways of evaluating the decay process, one by means of 
Feynman diagrammatic technique and the use of meson and nucleon propagators, the 
other one that uses operators and wave functions and we show the equivalence of the 
two methods. 

2.1.1 The propagator method 

The starting point is the A — > ttN Lagrangian, accounting for this weak process, 
which is given by 

C = -iGpii^^^iA^ + B^^,)r$^^jjA (?) + h.c. , (1) 

where Gf/x^ = 2.21 x 10^^ is the weak coupling constant. The empirical constants 
— 1.05 and B^^ — —7.15, adjusted to the observables of the free A decay, determine 

the strength of the parity violating and parity conserving amplitudes, respectively. 

The nucleon, A and pion fields are given by ■j/'n, V'a and 0^, respectively, while the 

isospin spurion (?) is included to enforce the empirical A/ = 1/2 rule, according to 

which the strength for the A — > 7T~p decay is double than that for A — > 7r°n. 

A practical way to evaluate the A width in nuclear matter and to introduce 

medium corrections is to start from the A self-energy, E, associated to the diagram 

of fig. 1, and then use the relationship 

r = -2ImS. (2) 
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Figure 1: Feynman graph for the free A self-energy of eq. (3). The A — > vrA^ "cut" is 
shown (dotted hne). 



Following standard Feynman rules, and making a non relativistic reduction of the 
operators in eq. (1), the self-energy is readily evaluated as 



- zm = 3{G,f^r J 4^.G{k - q)D{q) 



(3) 



where G and D are the nucleon and pion propagators, respectively, S = A^^, and 
P/fi = Bt^/2M with ju and M the pion and nucleon masses. 

The free nucleon and pion propagators are given respectively by 



G{k) 
D{q) 



1 



k^-E{k) + ie 
1 



(4) 



where the typical nonrelativistic approximations are done in the nucleon propagator: 
only the positive energy part of the nucleon propagator is taken and the factor M/ E 
{E being the relativistic energy) is set equal to unity. 

By using the free nucleon and pion propagators one obtains immediately the free 
A width |,|131 



free — J- A 



3{Gp^ 



2\2 



Sq 1 
(27r)3 2cj(g 



-2-K5{Ef,-uj{q) - E{k- q)) 



(5) 



where i?A, ^{q), E{k — q) are the energies of the A, pion and nucleon, respectively. 
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In a Fermi sea of nucleons, both the nucleon and pion propagators are changed to 



^/ N ^ 1 - n{p) nip) 

^^^^ " g0 2_^2_^2_n(g0^q) ' 

where Vat is the nucleon potential, n(g'', q = | g |) is the pion self-energy in the nuclear 
medium and n{p) is the occupation number in the Fermi sea, n{p) = 1 for \p\ < kp 
and n{p) = for \p \ > kp, with kp the Fermi momentum. An appropriate choice 
for Vn is the Thomas Fermi potential. Vat = —kp/2M, which becomes r dependent 
when the local Fermi momentum kpir) = (|7r^p(r))^/^ is used. 

The practical way to perform the integral in eq. (3) is through the Wick 
rotation displayed in fig. 2, where the analytical structure of the integrand is shown. 
The shaded region accounts for the discontinuity of the pion propagator due to Im 11. 
In one includes particle-hole (ph) and A-hole (Ah) excitation as a source of 
the p-wave pion self-energy plus and extra s-wave self-energy. Coulomb effects are 
neglected. Thus 

U{q) = U^^\q) + ffl{q) 
P 

-,U{q) 

with n(g) = -^^2 (8) 

1 - g'^uiq) 

and U{q) = UN{q) + U^{q) . 

In eq. (8), U{q) is the Lindhard function which gets strength from ph excitation, 
Up^{q), and Ah excitation, U/^{q), and g' is the Landau-Migdal parameter. In the 
A decay case, U^^q) has both real and imaginary parts, the latter one due to those 
contributions in the integration over the internal variables for which the ph excitation 
appears on the mass shell. The function Ui^{q) has also real and imaginary parts, the 
latter one coming from the (energy dependent) A width, which is relatively small in 
the A decay case since the pions emerge with little energy. Then, for practical purposes 
one can neglect Imf/^. In this case, the analytical structure of the integrand in eq. 
(3) with the medium propagators (6) and (7) shows, in addition to the cut (region of 
discontinuity in the real axis due to Imf/Ar 7^ 0), a pole in g° = Cj{q) corresponding 
to a renormalized pion energy for which 

u{qf - f - _ n(cD(g ), q) = . (9) 

Missing in fig. 2 is the pole of G{k — q ) corresponding to the second term in eq. (6). 
This pole lies in the lower half-plane of the figure and would contribute in the Wick 
rotation only when it happens to be in the third quadrant, i.e., k^ — E{k — q) — Vn < 0. 
But this corresponds to {k — q) very large where n{k — q) = ^ and hence this term 
does not contribute. 
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Figure 2: Analytical structure of the integrand of eq. (3) in the complex plane with 
the nucleon and pion propagators of eqs. (6), (7). The renormalized pion propagator 
pole uj{q) is shown. The dashed lines close to the real axis indicate the analytical cut 
from Imn(g°,g) related to the non-mesonic A decay channel. 



We come back to the evaluation of S(A;) in the nuclear medium, eq. (3), using 
the nucleon and pion propagators of eqs. (6) and (7). In the integral over g° along 
the contour of fig. 2, the contributions of the arcs at infinity vanish and the integral 
along the imaginary axis contributes only to Re E, given the property of the integrand 
= J*(g°*). Thus, only the residue of the integrand at the pole of the first 
quadrant is responsible for ImS and we obtain for the width ||5|,p^ 



QiGpfi')' J - n{k - q)]eik' - E{k - q) - Vm] 



(2vr) 



X 




Im 



g0 2 _ ^^2 _ ^2 _ n(gO, q) 



• (10) 



qO=kO-E{k-q)-VN 



In the discussion here we neglect the role of correlations and form factors, which 
are important only in the non-mesonic decay. We shall discuss this in detail when 
studying this decay mode in Sect. 3. 

In eq. (10) we observe the Pauli blocking factor, 1 — n. Since a A with momentum 
/c = decays into a nucleon and pion with q ~ 100 MeV/c, this momentum is 
smaller than the Fermi momentum for nuclear matter density, kp = 270 MeV/c, and 
the decay is forbidden by Pauli blocking, i.e, 1 — n{k — g) = 0. In finite nuclei it 
is still possible to have mesonic decay since the A wave function has some overlap 
with the nuclear surface where the Fermi momentum will be smaller than 100 MeV/c. 
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Moreover, the momentum distribution of the A leads to some spreading in the nucleon 
momenta allowing some of them to overcome the Pauli blocking. Nevertheless, the A 
mesonic width decreases drastically as a function of the mass number. 
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Figure 3: A self-energy diagrams included in eq. (3) with the nucleon and pion 
propagators of eqs. (6), (7). (a) Free self-energy graph, (b), (c) Insertion of p-wave 
pion self-energy at lowest order, (d) Generic RPA graph from the expansion of the 
pion propagator in powers of the pion self-energy, (e) s-wave pion self-energy at lowest 
order. The cuts represent the non-mesonic decay channel. 




A 



k-q ,1 



A 




(a) 



(b) 



Figure 4: Free and lowest order A self-energy graph. The dotted cuts represent the 
mesonic decay channel. 

The language of propagators used here is very adequate to provide a unified pic- 
ture of the A nuclear decay. Indeed, eq. (10) contains not only the modified mesonic 
channel but also the non-mesonic one. This can be seen diagrammatically by ex- 
panding the pion propagator and taking a ph and Ah excitation to account for the 
pion self-energy, 11. This is depicted in fig. 3. The imaginary part of a self-energy 
diagram is obtained when the set of intermediate states cut by a horizontal line are 
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placed simultaneously on shell in the intermediate integration. In fig. 3 we observe a 
source corresponding to placing on shell a nucleon and the ph of the pion self-energy. 
This corresponds to a channel where there are no pions and only nucleons in the final 
state. The physical process which has occurred is AA^ NN, which corresponds to 
the standard non-mesonic channel. Technically, it would be obtained by substituting 
in eq. (10) 

1 Im Uph 



Im 



gO 2 _ ^^2 _ ^2 



n 



2 



— /i^ — n| 



where Uph is the pion self-energy due to ph excitations. There is no overlap between 
Imnp/i(g°, q) and the pion pole in the propagator of eq. (10) and thus the separation 
of the mesonic and non-mesonic channels can be done in a clean way. 

The mesonic channel would correspond to a different cut, the one where the N 
and the tt are placed on shell. This is shown in fig. 4, where the diagram of fig. 4(a) 
corresponds to the lowest order in the expansion of the pion propagator. This term, 
together with the one in fig. 4(b) and further iterations contained in eq. (10), lead 
to a renormalization of the mesonic width, and an appreciable one, as it was shown 
in ref . |l^ . 



Technically, the mesonic width can be calculated either by subtracting the non- 
mesonic width from the total width, eq. (10), or, equivalently, by replacing the pion 
pole contribution in eq. (5) by the renormalized pion pole given in eq. (9). In this 
latter case the mesonic width is obtained, by analogy to eq. (5), from 

d^q 1 



3{Gpjj 



2\2 



X 




m 



2'k5{Ea - Cu{q) - E{k -q)-Vi 



N. 



(12) 



The qualitative reason for the drastic change on the mesonic width is given in 
|^,|l^]: the attractive character of the pion self-energy leads to a larger pion momentum 
for the same pion energy and thus, to a larger nucleon momentum by momentum 
conservation. Thus, the nucleon has more chances to have a momentum bigger than 
the Fermi momentum, therefore increasing the mesonic width. 



The width in finite nuclei is obtained in via the LDA 



r = y rfV|0A(f)pr(fc,p(f)) , (13) 

where (t>A is the A wave function. A further average over the momentum distribution 



of the A wave function is also done in |T3 . 



2.1.2 Finite nuclei approach. The wave function method 

The finite nuclei approach to the mesonic width was sketched in and carried out 
in detail in |jl5[|. The mesonic width is given, in analogy to eq. (5), by 
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(14) 



where (pN the wave function of the nucleon states and (p^^* corresponds to an outgo- 
ing solution of the Klein Gordon equation normalized to a plane wave asymptotically 
(e"*"^^). The index a stands now for tt^p or ir^n decay, with C^p^ = 4, C(") = 2, which 
one separates here since, due to shell effects, these channels can depart drastically 
from the elementary A J = 1/2 rule. 

The sum in eq. (14) runs only over non occupied nucleon states in the shell model. 
On the other hand, 0^"-** is a solution of the Klein Gordon equation with a proper 
optical potential (or pion self-energy, 11 = 2ujVopt), including Coulomb effects through 
Vc, i.e.. 



V' + /i' + 2uVoptix)] 0(-)(g, x)* = [iu~ V,{x)f <Pl\q, x)* . 



(15) 



The effects of using 0^^^* instead of a plane wave are rather drastic and increase 
the mesonic width in about two orders of magnitude in heavy nuclei [^,[1^,|T^, in 
qualitative agreement with the nuclear matter results of ref. [0. 

The arguments for the renormalization are expressed now in the alternative lan- 
guage as follows: the attraction caused by the pion self-energy leads to higher mo- 
mentum components in the pion wave function and, as a consequence, the matrix 
element of eq. (14), which involves the A wave function in the lsi/2 ground state and 
a nucleon wave function in an unoccupied orbit, is considerably enhanced. To see this 
let us write the pion wave function in terms of its Fourier transform 



ct>i-\lxy = J d'q'e-^'^ >.(g,g'), 

and recall that for | x |— > oo it behaves as e"*''^. For values of x inside the nucleus, 
due to the attractive character of the pion nucleus potential, the pion wave function 
builds up q' momentum components such that | q' \>\ q \. Then, the strength of 
the transition, | {N \ e"*'' ^' ^ 1 — iq'x — {q'xY /2 + . . . | A) |, appearing in eq. (14), 
is governed by the momentum q' which, being larger in magnitude than g, leads to 
enhanced matrix elements with respect to those using the asymptotic free pion wave. 
In the two languages the physical consequences are the same: an increased probability 
of reaching the unoccupied states and thus an enhancement of the mesonic width. 



2.1.3 Equivalence of the propagator and wave function methods 

Even if apparently the two methods discussed above look quite different, the physics 
contained in the propagator and wave function methods is equivalent. A formal 
derivation of the connection between the methods can be found in |TB|, which we 
reproduce here for completeness, complementing some points concerning the pion 
propagator in infinite nuclear matter and finite nuclei . 
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Let us start from the pion propagator in finite nuclei written in coordinate space 



(16) 



where (f)n{x) are the pion wave functions in the nucleus and e„ their corresponding 
energies. Ignoring pionic bound states, which do not play a role in our problem, we 
can identify the pionic wave functions by the asymptotic momentum q. Hence their 
energy is given by uj{q) = [q"^ + jJ^Y^"^- The sum over the index n is then replaced by 
an integral over q as given below 



(Pq (l)^{q,Xi)(l)l{q,X2) 



(27r)3 El-uiqf + iT] ' ^^^^ 

For simplicity in the derivation we shall take the parity violating part of the 
width (the one providing the largest contribution to the mcsonic decay) and will not 
distinguish between 7r° or tt" decay. Hence, from eq. (14) we obtain 



^F-l (27r)3 2a;(5) 



2T:6{EK-EN-uj{q)) 



X 



d?x<i)f^{x)<i)l{q,x)<j)*j^{x) 



(18) 



which can be rewritten as 



X 



/(pq 1 
7^^^^jpr27r5(EA -En- uj{q))(j)n{q, Xi)(f)l{q, fa) , (19) 



(27r)3 2u;(g) 
or, by virtue of eq. (17), as 

Ts = 3{GfI^YS^ J2 J c^'2;ici3x201(fi)0jv(fi)(-2)ImD,(fi, X2; E^^E^- En) 



N^F' 



X e{E^ - EN)<pA{x2)<i)*N{x2) 



(20) 



Now, in order to connect with eqs. (10) and (13) one makes a local density approx- 
imation. In the first step one evaluates F for a slab of infinite nuclear matter and 
in the second step one replaces the width in the infinite slab by an integral over the 
nuclear volume assuming slabs of matter in each d^r of the nucleus with local density 
p{r) and with a probability of finding the A particle given by |0a(^)P- This last step 
is implemented by means of eq. (13). Hence we should see how we reproduce now eq. 
(10) when we assume in eq. (20) a slab of infinite nuclear matter. For this purpose 
we have to replace for the nucleon sector 



E 

N^F 



V 



d^p 
(2^ 



:i - nip)) 



12 



En ^ E(p} - Vn , (21) 

and for the A wave function 

0A(f) ^e^'^^ . (22) 

Now, in the infinite slab of nuclear matter the pion propagator of eq. (17) is replaced 

by 

/^3g Qtq{xi-X2) 

where II{Et^, q) is the pion self-energy, which is a function of p. Note that for values 
of Xi, X2 far away from the nucleus, eqs. (17) and (23) are equivalent since there p = 
and n (in the LDA) will be zero. At other densities, 11 will be different from zero 
and the integral of eq. (23) gives rise to other momentum components, modulating 
the plane wave in the numerator and providing a kind of WKB approximation to the 
wave functions in the numerator of eq. (17). The LDA gives rise to a variable local 
momentum and hence a distorted pion wave. 

By substituting eqs. (21), (22), (23) in eq. (20) we obtain: 



r, = -6(G^^^)25^|-g_y-^(l-n(p-))lmZ^.(g)%°) 



qO=EA-E{k-q)-VN 



X J d^xi£x2^e'P^^'-^^h'^^"^^e-'''^"^^ , (24) 



with Dj,{q) given by cq. (7). Finally, by means of the relationship (27r)^5^(0) 
/ d^x — V we can cast eq. (24) as 



r5 = -eiGpp'fs' 



d^q 
(2^ 



n{k-q))lmD„{q)e{q^) 



, (25) 



qO=EA-E{k-q)-VN 



which coincides with the s-wave contribution to F from eq. (10). This establishes 
the equivalence between the finite nuclei wave functions method and the propagator 
method, the latter one complemented with the LDA. 

A further consideration can be done regarding the pion propagators in finite nuclei 
and nuclear matter, eqs. (17) and (23) respectively. Neglecting, for simplicity, the 
Coulomb potential, assuming a smooth energy dependence of the pion self-energy and 
taking 11 for a value of E^^ around the pion pole, the Klein-Gordon equation leads to 



-V^ + p' 



+ Uiq,p{x2))-El 

d^q uj{qY-El 
(27r)3 El - uj{qf 



D^''{xi,X2,E^) = 
^{q, xi)(f)l{q, X2) = -S^{xi - X2) 



(26) 
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where in the last step we use the fact that the pionic wave functions form a complete 
set of states. 

We can proceed equally with the nuclear matter propagator and we obtain 



-V2 ^ ^^ + T{{q,p{x2))-El\D^^'{xi,X2.E, 



i2nrEl-q^-p^-n{q,p{x2)) 



ig{s,-x,) ^ _^3(^^ _ _ (27) 



The integral over q gives the same results in eqs. (26) and (27), leading to a local 
function, and providing further insight in the meaning of the LDA approximation in 
inclusive processes. 

Note also that for a fixed E^^ of the pion, the value of q at the pole of the two 
propagators is different, but also its meaning. In the case of finite nuclei we have 
Et^ = uj{q) and the value of g is the asymptotic one when the pion leaves the nucleus. 
In the case of infinite nuclear matter we have E'^ = u;(g )^ + n(£^7r; p(^)) which gives 

En ^ uj{q) + = uj{q) + Kpt(r) . (28) 

Hence, in this latter case, the value of q is the local momentum of the pion such that 
the total energy, kinetic plus potential, is the asymptotic energy of the pion. 

2.2 The mesonic width and the occupation number 

We have seen that Pauli blocking is the major factor in reducing the mesonic width 
of heavy A hypernuclei. It was suggested that because real interacting nuclei have 
the "occupied" states partly unoccupied, the mesonic width should be enhanced with 



respect to a calculation with fully occupied states up to the Fermi level [0. In the 
nuclear matter approach of Sect. 2.1.1 this is easily visualized by recalling a realistic 
picture of the occupation number of the Fermi sea [^, which is depicted in fig. 5. 



For the states below the Fermi momentum the level of occupancy is of the order of 
85% and, by assuming that in the A decay the nucleons can occupy the 15% vacancy 
of these states, we would guess that the mesonic width would stabilize at the level of 
about 10% of the free width for heavy nuclei (taking into account the absorption of 
pions on their way out of the nucleus). If this were the case the mesonic width could 
serve as a measure of the occupation number in the Fermi sea. The argument is very 
appealing and intuitive, however, it is incorrect and leads to an overestimate of the 
width by about three orders of magnitude in heavy nuclei. 



The detailed discussion of this problem was done in ref. [^. The falacy in the 
argumentation lies in the the following: since the Fermi distribution n{p) of eq. (6) 
is only a first approximation to the momentum distribution, it looks like a sensible 
improvement to substitute the nucleon Fermi propagator of eq. (6) by 

l-ni{k) ^ niik) ^^^^ 



kO-E{k)+ie k^ - E{k) - ie 
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n(k) 



Figure 5: Schematic representation of the nucleon occupation number for an inter- 
acting Fermi sea. 

where ni{k) is the reahstic occupation number in nuclear matter and is then ^ 

— * 

for all values of k. This is actually and approximation often found in the literature, 
particularly in works of condensed matter and heavy ions. Given the analytical 
structure of the propagator, eq. (29), identical to the one of eq. (6), the formal result 

— * — * 

for the width is the same as in eq. (10) replacing only 1 — n(k — g) by 1 — ni{k — g ), 
in which case we would find widths of the order of 15% of the free width. However, 
eq. (29) is not an improvement over the propagator of eq. (6), but quite the opposite, 
as we shall see. The realistic nucleon propagator for an interacting Fermi sea is given 
in terms of the spectral functions by 

G(*;», 4) = r du,-iP^ + r <i.#<^^ , (30) 

^ ^ i-oo k^-u-ie J Hp k^ -u + ie ^ ' 

with /ip the chemical potential. 

When performing the calculations of the mesonic width with this nucleon propa- 
gator one obtains the factor 

/ duSp{u,k- q)2T:5{k^ -u -u{q)) (31) 
J tip 

replacing the factor 

[1 - n{k - q)]2T:5{k^ - E{k - q) - uj{q)) (32) 

in eq. (5), when the Pauli blocking factor of eq. (10) is implemented. Eqs. (31) and 
(32) bare some intuitive resemblance because 

I du!Sp{u!, k — q) = 1 — / Sh{uJ, k — q)duj = 1 — ni{k — q) . (33) 

J l^p J — oo 



I 

I 
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However, the presence of the 5 function in eq. (31) prevents from factoring out the 
integral of 5*^ shown in eq. (33). Furthermore, because of restrictions of the phase 
space (energy and momentum conservation) the range of values of u allowed are very 
small compared to the range (/ip, oo) needed in eq. (33) to obtain 1 — ni{k — q) of the 
interacting Fermi sea. In physical terms we can interprete it in the following way: the 
occupation number ni{k — q) is an integral for all the energies of the nucleon, u, of 
the probability of finding a nucleon with momentum k — q and an energy u, which is 
given by the spectral function Sh{uj, k — q). However, in a physical decay process we 
have conservation of energy and momentum and hence there are severe restrictions 
to the values of the energies that the nucleon can have. This is why the occupation 
number riiik — q) cannot be factored out. 

The actual calculations carried out in ref. |]2T]] showed that for light and medium 
nuclei the use of the spectral representation for the nucleon propagator, eq. (30), 
instead of the one of the noninteracting Fermi sea, eq. (6), has negligible consequences 
in the mesonic width (of the order of 6% corrections in O). The corrections can be 
of the order of 50% in heavy nuclei, but in all cases, when the pionic renormalization 
is taken into account, one can disregard these effects. 

These findings have been of relevance in showing similar problems in the study 
of other physical processes, like in the contribution of the pion cloud to nucleus 



scattering where it is shown [22| that one cannot relate the effect of the pion cloud to 



the pion excess number in the nucleus as assumed in refs. Similar spectacular 

differences between the use of the occupation number and the spectral function are 
seen in the study of the nuclear structure functions of deep inelastic scattering at 
values of X > 1 |2^. For values of x = 1.3 the use of the occupation number gives 
values of the structure function two orders of magnitude bigger than the results with 
the spectral functions, which agree with experiment [p^] . 

2.3 Results for the mesonic width 

We separate here three regions of heavy, medium and light nuclei, where different 
physical phenomena are explored by means of the mesonic decay. 

Experimental data for light and medium nuclei of the last decade can be found 
in 0,|2|,E9|,Rg,pl. 



2.3.1 Medium and heavy nuclei 

In refs. IIT5|,|IB101 one can find abundant results in different nuclei which are rather re- 



alistic. These results have been recently improved ||3^ by a more accurate description 



of the energy balance in the particular reactions, taking into account transitions to the 
bound and continuum nucleon states and using a pion nucleus optical potential which 
has been derived theoretically and leads to a good description of pionic atoms data 
and of elastic, reaction and absorption cross sections in the scattering processes [p3| . 
The imaginary part of the pion nucleus optical potential is split into two terms related 
to pion absorption and quasielastic scattering. In P2| the pion quasielastic events are 



16 



not removed from the pion flux, as it corresponds to the actual experimental observa- 
tion, while the use of a full distortion of the pion with the total optical potential, as 
done in ||l^,|^,|l3, inevitably removes the pion quasielastic events, together with the 
pion absorption events. Though conceptually important, this refinement turns out 
to be of little practical relevance in the present problem given the small energy that 
the pions carry and the very small phase space for quasielastic collisions P^. How- 
ever, other considerations, particularly the energy balance in the reactions, makes the 
widths in heavy nuclei for 7r~-decay about one order of magnitude smaller than those 
of ref. ilT 
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Figure 6: Pionic decay rate for 7r° and vr^ as a function of the mass number (of 
the host nucleus, ^^O, ^°Ca, ^^Zt, ^^^Ba, and ^°^Pb). The two lower lines show the 
calculation with plane waves for the pion and the two upper lines the results with 
pion distorted waves. 



In fig. 6 we show the prediction of ref. [32| for different nuclei and for 7r° and tt^ 
decay, with plane waves and the renormalized pion wave function. The drastic effects 
of the pion renormalization are seen there and are a bit smaller than in former works 
because the energy balance forces the pions to come out with smaller energies than 
in the previous approaches and the attractive effects of the p-wave part of the optical 
potential are then diminished. Shell effects and energy balance are important here 
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because of the small energies left for the pions. They are responsible for the slight 
increase of the A width in j^^^Ba with respect to ^Zr seen in the figure for the plane 
wave calculation. When the interaction of the pions with the nucleus is considered, 
the extra attraction felt by the pions gives more room for pionic decay such that the 
nuclear energy differences in the energy balance are not so crucial. 



2.3.2 Shell effects in medium nuclei 

The energy balance for some nuclei makes the results somewhat surprising. Of par- 
ticular relevance are the results in shown in Table |l[ 



Table 1: Pionic partial decay rates in the decay of 





F^o/Fa 


r.-/FA 


F^o/F^- 




0.159 
0.13 


0.086 
0.098 


1.86 
1.32 


EXP: [pop 
EXP: [31 


0.217 ±0.084 


U.UOZ_o.035 





In ref. [|l^, using different pion optical potentials, the authors quote 0.169 and 
0.126 for F^o/Fa and 0.134 and 0.107 for F^-/Fa. Although with large errors, the 
experimental results confirm these striking theoretical predictions which show a large 
departure from the A/ = 1/2 rule in nuclei. Indeed, we would naively expect 
Ftto/Ftt- ~ 0.5 (this should be the case for N = Z core nuclei), while it is of the 
order of 2 for the results of ref. [^ and slightly larger than one for the results of 
refs. []16|,0- This is mostly due to nuclear shell effects because of the large differ- 



ences between the Q values met in these reactions. 



2.3.3 Light nuclei: Short range repulsion and quark models 



Another interesting finding is seen in very light nuclei. The mesonic width of \B.e 
has attracted particular attention. There, in addition to the pion renormalization, 
the repulsive character of the AA^ interaction and the relatively weaker medium range 
attraction, compared to the A^A^ interaction, push the A wave function to the surface 
of the nucleus, weakening the Pauli blocking effect and thus enhancing the mesonic 
decay [3i,35|. The experimental numbers clearly favour potentials with a repulsive 
AA^ core. One should note that such a repulsion automatically appears in quark based 



models of the AA^ interaction [Q . The study of the \}ie decay using a quark model 



based hypernuclear wave function leads to the results shown in Table y, where 
they are compared to those of [Q, which uses a A wave function from the modified 
YNG AA^ interaction [39| with a strong short range repulsion, and to those of [HH] for 



the Isle A-nucleus repulsive potential at short distances [^ . 

There is a fair agreement of the theoretical results []37| , p8| , ^] with experiment |pO|| , 
account taken of the large experimental errors. 

Although the effect of the renormalization of the pion in the medium leads to 
an enhancement of the inclusive mesonic rate, it is possible to find particular decay 
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Table 2: Pionic partial decay rates in the decay of ^He 









F./Fa 




m 


0.431 


0.239 


0.670 


m 


0.393 


0.215 


0.547 




4C] 


0.40 


0.20 


0.5 


EXP: [0 


0.44 ±0.11 


0.18 ±0.20 


o.59i^:M 



channels where the effect is reversed. This has been shown in ref. for the decay 
aH —y^ He ± 7r~. The final nucleus ^He selects the s-wave part of the pion optical 
potential, which is repulsive, hence leading to a reduction of the rate with respect to 
a calculation using free pion waves. This also means that, through the study of the 
mesonic decay in selected channels, one can learn about different parts of the pion 
nucleus optical potential. 



The mesonic decay of the hypertriton, aH, has also been studied recently 



using wave functions for the hypertriton, 3N bound and 3A^ scattering states which 
are solutions of the Faddeev equations. The ^He±7r~, p+d+TT~ and p+p+n+TT~ decay 
channels and the corresponding ones for 7r° have been studied. The total mesonic 
width is 92% of the free A width and compares favourably with the experimental 
data for the hypertriton lifetime, which ranges between (2.20iJ:i) X 10 sec to 
(2.64l°-^I) X 10-10 sec fpl. 



-0.54; 

An interesting problem related with the decay of light hypernuclei is the ratio 
tt^/tt- in the decay of ^He. This ratio is measured to be about 5% and cannot come 
from the direct decay of the A (A — >■ 7r~p, ir^n). The problem concentrated some 
attention in the past where rescattering of the vr^, vr*^ + p — >■ vr"*" + n, or E"*" 
decay, S+ n^n, following the Ap S+n conversion were considered. However, 
the results gave a fraction much too small compared to experiment. The problem 
was revived recently using updated information, and a result for that fraction of 
1.2% was obtained, which, although being a factor of two larger than in [^], is still 
too small compared with the experimental value. 

A more recent work [^] proposes a solution which consists of adding to the AHe 
wave function some S"*" ± ^H component, which is relatively sizable in view of the 
strong AA^ J^N conversion. The S"*" would then decay into H'^N —>■ TT~^nN, which 
leads to a s-wave spectrum as observed experimentally 

On the other hand, a different view is taken in , where it is suggested that the 
A/ = 3/2 component of the weak interaction is the responsible for the relatively high 
fraction of vr"*" emission. 



2.4 Conclusions on mesonic A decay 

We have made a review of the present situation concerning the mesonic decay of A 
hypernuclei. We have established the formal link between the propagator method, 
where the huge enhancement of the pionic decay width was first reported, and the 
finite nuclei approach with wave functions and matrix elements. Shell effects and 
precise values of the nuclear binding energies are also important in the mesonic width 
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and they are best taken into account in the finite nuclei approach. The intuitive and 
appeahng, but falacious, hnk between the nucleon occupation number and the mesonic 
width has also been discussed, which has served to unveil rough approximations used 
in other processes to link the pion excess number with contributions of the nuclear 
meson cloud to some physical observables, hke K'^ nucleus scattering or deep inelastic 
scattering. We have also discussed the relevance of the short range KN repulsion in 
the mesonic width of light hypernuclei and showed how different models, all of them 
accounting for this repulsion, can provide a fair description of the experimental data. 
We should note in this respect that microscopical quark models lead naturally to such 
kind of repulsion, although a description in terms of the more conventional meson 
exchange picture is also posible. 

With the limited amount of experimental data available on the mesonic channel, 
the amount of physical information obtained is remarkable. As just mentioned, there 
is support for strong KN repulsion at short distances. The process also provides us 
with the most striking renormalization effect due to the pion nucleus interaction. The 
sensitivity of the A decay to the pion nucleus optical potential can serve as a tool to 
choose between different theoretical descriptions of the complex mechanisms of pion 
nucleus interaction. The decay channel into tt'' can be an excellent instrument to 
learn about tt*^ nucleus interaction, and so on. 

It is clear that a systematic experimental search in many nuclei of the mesonic de- 
cay channel will provide us with very valuable information to unravel the intricancies 
of the pion nucleus interaction or the elementary properties of the KN interaction, as 
well as proper nuclear structure details of the A hypernuclei themselves. 
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3 Weak non-mesonic decay of A hypernuclei 



The pionic decay mode of A hypernuclei is neghgibly small in medium to heavy 
hypernuclei due to Pauli blocking acting on the final nucleon which emerges with a 
very small momentum. However, the medium effects are responsible for the opening 
of new decay channels, where there are no pions in the final state. These non-mesonic 
decay modes of A hypernuclei can be viewed as coming from the mesonic decay when 
the pion is produced in a virtual state and absorbed by one or more nucleons. The 
one-nucleon induced decay mode corresponds to the AA^ NN transition and can be 
interpreted as being mediated by the exchange of one pion (or more massive mesons). 
In the two-nucleon induced decay of nuclei, the transition ANN NNN can be 
interpreted as coming from the absorption of the virtual pion by a pair of nucleons 
that are correlated by the strong interaction. No matter what the interpretation of 
the different non-mesonic decay mechanisms is, what happens is that the mass excess 
in the initial state (176 MeV) is now converted exclusively into kinetic energy of 
the final nucleons, which may emerge with large momentum values not restricted by 
Pauli blocking. Hence, except for the very lightest ones, A hypernuclei decay mainly 
through the non-mesonic mechanisms. 



3.1 One-nucleon induced decay 

The decay of hypernuclei through the one-body induced mechanism AA^ NN offers 
the best opportunity to study the AS = — 1 nonleptonic weak interaction between 
hadrons due to the practical impossibility of having stable A beams. Since the typical 
momentum exchanged is around 400 MeV/c, the process is short ranged and will not 
be too sensitive to the nuclear structure details. The AA^ —>■ NN reaction, similar 
to the AS* = A^A^ — > A^A^ weak reaction but in the AS* = —1 sector, allows 
one to study not only the parity violating (PV) part of the interaction, but also 
the parity conserving (PC) one which, in the A^A^ — A^A^ case, is masked by the 
strong interaction. Note, however, that the study of the inverse reaction pn pA 
with moderate energy and high intensity proton beams could be feasible at RCNP 
(Osaka) g or COSY (Jiilich) 0. 

The non-mesonic decay of hypernuclei has been studied for more than thirty years. 
The experimental and theoretical status on the one-nucleon induced hypernuclear 
decay was extensively reviewed by Cohen in 1990 [Q. For this reason we will just 
update the present situation in the next paragraph and only the new achievements 
made since 1990 and how they compare with previous works will be covered in more 
detail in the coming sections. 

Block and Dalitz built a phenomenological model to describe the decay of 
s-shell hypernuclei in terms of a few elementary spin-isospin reaction rates. Although 
no assumptions were made on the dynamics of the AA^ NN reaction, this model 
has been (and is) extremely useful because, by fitting to the empirical data for light 
hypernuclei, one can extract information on the characteristics of the reaction, such 
as the validity of the A J = 1/2 rule. The first microscopic calculation was that of 
Adams who considered the one pion exchange mechanism to describe the non- 



21 



mesonic decay of a A in nuclear matter and found a large sensitivity to the short 
range correlations (SRC) induced by the strong interaction in the initial AA^ and 
final A^A^ states. A unified treatment of the mesonic and non-mesonic channels based 
on Feynman diagrams and propagators was developed by Oset and Salcedo [|13[, where 



the medium effects in the renormalization of the pion were considered. Being a short 
range process, it was clear that the one pion exchange (OPE) mechanism might be 
insufficient to explain the data. McKellar and Gibson ||5^ made the first attempt to 
include the heavier mesons by developing a model for the exchange of the p meson. 
This was applied later to the decay of light hypernuclei |^ but the results depended 
on the relative sign between the tt and p mechanisms that was not fixed by the model. 
Based on a pole model for the PC vertices and on soft meson theorem techniques for 
the PV ones, Dubach et al. constructed a full one meson exchange potential and 
preliminary results for nuclear matter were obtained. The details of the model and 
final results have recently become available |^ . An extension of the model including 
form factors and some exchange diagrams related to the antisymmetrization of the 
final two-nucleon state has been recently applied to the decay of finite hypernuclei 
5^. The exchange of a and p mesons has also been considered within the image 
of correlated two pion exchange, where A and S are excited in the intermediate 
state [^,^,|6l],^]. A description of the AA^ — A^A^ decay in terms of quark degrees 
of freedom was first attempted in refs. p3[, where a hybrid model that combined a 



quark Hamiltonian at short distances and OPE at long distances was employed. Oka 
and collaborators combined a direct quark mechanism with a one pion exchange 
contribution for the description of light hypernuclei. A more realistic A wave function 
was used in a more recent work where also the phase between the direct quark 
and the pion mechanisms was established. One kaon exchange amplitudes, which 
partially cancelled the pion ones, were also calculated in ref. |^B[ together with the 
quark contributions. All the models and results of the recent works will be extensively 
reviewed below. 



3.1.1 OPE results 

The OPE mechanism for the AA^ NN transition displayed in fig. |^ involves a 
weak and a strong vertex. 

As already shown in eq. (1), the weak Hamiltonian is parametrized in the form 

TiZ^ = iGfp''^^{A^ + 5^75)f0,z^A (0) , (34) 

which implements the A/ = 1/2 rule through the isospin spurion (5). The con- 
sequences of breaking this rule in the non-mesonic decay of hypernuclei have been 
recently explored and will be presented in Sect. 3.1.4. 
For the strong vertex the pseudoscalar coupling 

^NNTT = ^5'NN7rV'N75^07rV^N , (35) 

has usually been taken since most calculations are nonrelativistic and this vertex is 
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Figure 7: One pion exchange diagram for the AA^ NN transition 
equivalent to the pseudovector one. The nonrelativistic transition potential reads 

= -G./i^l^ (a. + f=a,g) ^P^.^ (36) 

where q is the momentum carried by the pion directed towards the strong vertex, /i 
the pion mass, M the nucleon mass and M the average between the nucleon and A 
masses. 

By introducing the tensor operator S'i2(r) = 3(aig)(a2g) — (?i(T2, the OPE transi- 
tion potential can be splitted into central, tensor (adding to a total PC part) and PV 
pieces. The decay rate can be obtained in terms of < / | 14- I ^ >, the expectation 
value of the potential between initial AA^ and final A^A^ two-body states. 

There are a few relativistic corrections which can be implemented. The vertices 



can be made relativistic (see eq. (32) of [^] and eq. (3) of ||68|) and the pion 
propagator non-static, i.e., — >■ — g^. Furthermore, one should use, for consistency, 
relativistic wave functions |]69|, [70|] . Both types of corrections are of the order of 10% 



and of opposite sign, so we shall ignore the relativistic effects in what follows, 
a) Nuclear matter 

The first work to consider the OPE mechanism was that of Adams ^3| , where the 
importance of strong and final state interactions in the AA^ NN transition was 
first pointed out. The calculation was made in nuclear matter and only the L = 
AN relative motion was retained. The same approach was followed by McKellar 



and Gibson who made the first attempt to improve on the OPE mechanism by 
including the p meson as well. This will be discussed in the next section. Their 
uncorrelated one-pion rate was consistent with that of Adams, after correcting the 
too small coupling constant, but the effect of correlations was found less strong. This 
was confirmed by the two other nuclear matter calculations that appeared at that 



time. Oset and Salcedo [Ol used a Green's functions formalism to treat the mesonic 



and non-mesonic channels in a unified way. Their method includes automatically all 
the partial waves of the relative AA^ motion. The OPE mechanism was also studied 
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by Dubach and collaborators in the context of a one meson exchange model for the 
AA^ NN transition 

The results for the one-nucleon induced non-mesonic decay rate in nuclear matter 
are summarized in Table |^. The values of this table, and all others referring to 
decay rates, will be given in units of the free A width. Fa = 3.80 x 10^ s~^. The 
uncorrelated results obtained by the various groups are very similar. Substantial 
diferences are found when SRC are considered, although Adams used an unusually 
strong tensor correlation that, when omitted, would bring the correlated result up to 



1.57. The discrepancy of a factor of 2 between refs. |]13[ and when form factors 
and short range effects tied to the initial and final strong correlations are included, 
can be understood from differences in the model ingredients. In ref. a combined 
form factor of the type $(g) = (A^ - yU^)/(A^ + g^) is employed with A^ = 20yU^, 
while ref. |TB] uses this monopole form factor at each vertex with A^r = 1300 MeV. 
For the relevant momentum transfer of 400 MeV/c, the ratio of form factors squared 
is 1.8, while the different correlation functions account for the remaining difference. 
McKellar and Gibson use a Gaussian type, 1 — exp(— with r = 0.75, while Oset 
and Salcedo use a Bessel type, 1 — joiqc^), with = 3.93 fm~^. As seen in ref. 



the rates calculated with these two choices differ by about 30% and the Bessel type is 
closer to what is obtained from microscopic G-matrix calculations using the Nijmegen 
AA^ interactions. 



Table 3: One pion exchange contribution to the A non-mesonic decay rate in nuclear 
matter. 



Group 


uncorr. 


corr. (no FF) 


corr. + FF 


Adams 153 






3.47* 


0.38 




McKellar-Gibson 


54 




4.13 


2.31 


1.06 


Oset-Salcedo |13 




4.3 




2.2 


Dubach et al. 57 




4.66 


1.85 





* Corrected by a factor of 6.81 as pointed out in ref. |54|[ 



In summary, the uncorrelated OPE contribution to the A decay rate in nuclear 
matter turns out to be around 4 times the free A width and the results are lowered by 
a factor of 2 when appropriate short range correlations in the initial and final state 
are considered. 

b) Finite nuclei 



The existence of new experimental data for the non-mesonic decay rates p0|,|72 
and the possibility of obtaining polarization observables, such as the asymmetry in 
the distribution of protons emerging from the weak decay of ^^C []73[ motivated a new 



calculation of the non-mesonic decay performed directly for the finite hypernuclear 
system [^,[7^]. The formalism was chosen relativistic in view of the success of this 
model in reproducing polarization observables in proton nucleus scattering reactions. 
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The AA^ — > A^A^ transition was obtained from the relativistic Feynman amphtude 



tfi = iGpl^'g^^nJ c/W?/%^(x)(A+S75)^«Jx)A,(x-i/)/(| f-f |)%^(?/)75^'«^(x) , 

(37) 

where the initial A and nucleon were solutions of the Dirac equation with scalar and 
vector potentials, chosen to reproduce the nucleus form factor and the nucleon and A 
separation energies. The two emerging nucleons felt the effect of a relativistic optical 
potential. The short range effects were considered via a phenomenological correlation 
function, /, modifying the pion propagator A^^. The value of the free rate for 
was found to be about a factor of 2 smaller than the nuclear matter results. This is 
reasonable because a finite hypernucleus has low density regions which give a smaller 
contribution to the decay rate. It was surprising, however, that the rate was reduced 
by a factor of four due to SRC and FF at the vertices, an effect twice as large as that 
observed in the nuclear matter calculations. To adscribe the differences to the finite 
size of the system is not reasonable because correlations are short ranged and they 
should affect similarly the infinite and finite systems. Another possibility would be to 
relate the discrepacies to the relativistic model but ref. already showed that, in 
the absence of correlations, the relativistic and nonrelativistic treatments were giving 
the same results to within 10%. 
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Figure 8: Schematic model for correlations in the OPE potential. Bare OPE diagram 
(a) and OPE with simultaneous exchange of a heavy meson (b). 



The answer to this problem was partially given in ref. [|73[, where it was shown, 
through a nonrelativistic reduction of the correlated Feynman amplitude of eq. (|37|) , 
that one did not obtain the same correlated potential as the one used in the nonrela- 
tivistic approaches, namely /(/')Kr(^)- The final answer was given in ref. |]7l| where, 
through a microscopic model for correlations, it was shown that the correct proce- 
dure is that followed by the nonrelativistic models. The essential point is that the 
origin of correlations is tied to the exchange of a heavy meson, as e.g. the u meson, 
represented by the wavy line in fig. ^ Then the vertices involved in the weak transi- 
tion mediated by the pion (dashed line) connect final state spinors with intermediate 
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spinors which depend on the momentum over which one must integrate. This can be 
clearly seen from the momentum space expression corresponding to the PC piece of 
the weak transition potential 

D{q)u{ps)-f5u{p3 + q)u{p4h5u{p4 -?)••• (38) 



(27r) 



where D{q) is the pion propagator and the dots simbolize additional elements of 
the amplitude which are not relevant for the discussion. However, the relativistic 
procedure followed in refs. [^9|,^ was equivalent to working with external momentum 



spinors which, therefore, factorize out of the integral 



d'^q 

u{pz)iMvi)u{'Pa)i^u{p2) I ^^^^(^) • • • (39) 



In coordinate space, this procedure leads to oscillating potentials |7T|] which produce 
an artefactual reduction of the rates. 

The calculated non-mesonic decay rates of several hypernuclei have been reported 
recently in the framework of a nuclear shell-model and a nonrelativistic transition 
operator. Details on how the hypernuclear transition amplitude can be decomposed 



in terms of two-body amplitudes for nucleons in any shell can be found in [75|. The 
nucleons and the A are described with harmonic oscillator wave functions, SRC are 
included via a phenomenological correlation function inspired in G-matrix calcula- 
tions with realistic AA^ interactions and final state interactions (FSI) are treated in 
terms of two-body A^A^ scattering states obtained from realistic A^A^ interactions. A 
phenomenological form factor is included at each vertex as well. The effect of the 
different ingredients is visualized in Table | for the decay rate of ]^C, splitted into 
central, tensor (adding to a total PC rate) and PV contributions. 



Table 4: OPE contribution to the non-mesonic decay rate of from ref. [| 





FREE 


SRC 


SRC+FF 


SRC+FF+FSI 










phenom. 


Nijm93 


Central rate 


0.282 


3.4 X 10^^ 


1.3 X 10-^ 


4.2 X 10-3 


3.3 X 10-3 


Tensor rate 


0.858 


0.781 


0.637 


0.685 


0.566 


PC 


1.140 


0.785 


0.650 


0.689 


0.531 


PV 


0.542 


0.447 


0.389 


0.421 


0.353 


P/Pa 


1.682 


1.232 


1.038 


1.110 


0.885 



The free central rate is reduced drastically by SRC, however, most of its contribu- 
tion comes from the delta function which is completely eliminated by the correlation 
function, which is zero at the origin. The tensor rate is reduced by 10% by SRC 
and by 20-35% once FF and FSI are included. The PV rate represents 40% of the 
total rate, differing from previous nuclear matter results, which reported either negli- 
gible [54 1 or 15% 1 56 1 contribution to the rate. However, a simple estimate of the PV 
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contribution, relative to the PC one, for the most relevant momentum transfer (400 
MeV/c) gives 



2M 



0.57 , 



which yields a 35% fraction of the total for the PV rate. The last two columns 
reflect the importance of using a realistic A^A^ wave function in the final state. A 
phenomenological NN correlation function of the type 1 — jo{qcr), with Qc = 3.93 
fm~^, gives a rate 25% larger than that obtained from a NN wavefunction based on 
the Nijmegen93 potential The NN interaction is relevant because it modulates 
the wave function at the short relative distances and prevents the contribution from 
the short part of the interaction. One may wonder about the interaction of the 
emitted nucleons with the mean potential of the nucleus. This can be taken into 
account by means of a complex optical potential, the real part of which is moderate 
and has practically no effect in the transition rate. The imaginary part would induce 
a loss of flux due, in the present case, to A^A^ collisions. However, the net effect 
is a redistribution of the original A^A^ emission strength into different multinucleon 
channels |77| without modifying the original decay rate, which is the magnitude in 
which we are interested. 

One of the interest in doing finite nuclei calculations was to establish the role of 
relative L ^ motion in the non-mesonic weak decay. It was especulated that the 
p-shell rate would be basically coming from relative AA^ motion in P wave, due to 
supression of S-wave by the corresponding L = 1 center-of-mass wave function which 
is zero at the most probable back-to-back kinematics. However, it was found ||7^ that, 
after integrating over all possible angular possibilities, 90% of the p-shell contribution 
comes from the relative AA^ L = motion, which has been confirmed by the most 
elaborated model of ref. ||58[] . 

In Table § finite nuclei results from different works are compared. The results 
within brackets, if available, correspond to the calculation omitting SRC, FF and 
FSI, which will be referred to as the uncorrelated ones. 



Table 5: OPE contribution to the non-mesonic decay rate of hypernuclei. Results 
within brackets correspond to the calculation omitting SRC, FF and FSI. 





iHe 


A ^ 


Takeuchi et al. 
Dubach et al. | 

Parreno et al. 
Oset-Salcedo ; 


57 
58 

a| 


1 


(0.52) 0.14 
(0.6) 0.9 
(1.6) 

(0.98) 0.5 
1.15 


(3.4) 0.5 

2.0 |79| 
(1.68) 0.89 
1.45 


EXP: 


0.41 ±0.14 g 


1.14 ±0.2 Q 
0.89 ±0.18 |72| 



Compared to the other predictions, a quite low rate for ^He is found in ref. 



55|. 



Part of it is due to the use of a more realistic A wave function pushed to the surface 
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by the effect of the repulsive A nucleus interaction in light hypernuclei. But most 
of the reduced rate comes from a much stronger effect of FF and SRC as they use 



the model of ref. |^ discussed in the nuclear matter section. The recent results 
of Dubach et al. are somewhat confusing and the lack of computational details 
makes the comparison with other works difficult. However, some general points are 
worth mentioning. Their uncorrelated rate for j^^C is a factor 2 larger than that of |^ 
and the correlated one a factor 2 smaller, so the effect of correlations is substantially 
larger than that found in their nuclear matter results. It is also surprising that their 
correlated rate for |He is larger than the uncorrelated one and a factor of 2 larger 
than the result. An explanation could lie in the possibility of some misprints. 
Their rate quoted more than 10 years ago in some conference proceedings ||79| 



turns out to be around 2. Moreover, the preprint of their most recent paper p7 



quotes 1.6 for the uncorrelated \B.e result. This two corrections would lead to more 



consistent results and would differ only in a factor of two from those of ref. ||5^. We 



note that the model of ref. [|57|1 does not contain form factors at the vertices, which 



would provide some additional reduction. The LDA results of ref. [|T3| appear to be a 
little bit too high when compared to the direct calculations on finite nuclei. However, 
we have checked that, when using the same A wave function as in ref. (slightly 
more extended), a similar correlation function — controlled in |13| by the Landau g'j^ 
parameter — , and neglecting the pion renormalization, the LDA result lowers down 
to 1.05, nicely agreeing with the value 1.038 quoted in Table § without FSI effects, 
which are absent in |T^. This consistency check reinforces the idea that the LDA is 
an excellent tool to study processes that are not too sensitive to nuclear structure 
details or low energy nuclear excitations as is the case of the non-mesonic weak decay 
of hypernuclei. 

After this comparison, we can safely claim that, in general, the OPE mechanism 
predicts non-mesonic decay rates of finite nuclei in agreement with the experimental 
values. 

However, the OPE model is unable to explain the ratio r^/Fp between the neutron- 
induced rate [An — >■ nn) and the proton- induced one (Ap — > np) . In Table |^ several 
predictions for this ratio are compared with the experimental results. While the data, 
with large error bars, seem to suggest a value of around 1 or larger (except for the 
very light hypernucleus y^^He |^,|ST|), all models give values around 0.05-0.2. This is a 
well known characteristic of the OPE mechanism, which favors the tensor transitions 
over the central and parity violating ones. Since the AA^ pair is mainly moving in 
L = states and the tensor operator induces ^S'i(AA^) — ^Di{NN) transitions, 
the antisymmetric final system will be predominantly in an isospin 1=0 state, which 
only np pairs can achieve. Hence, the transition An nn is highly suppressed in 
the OPE model. These arguments have been more quantitatively stated in ref. P|, 
within the propagator method. In the limit of zero momentum for the A and small 
Fermi momentum for the nucleons, the ratio F„/Fp is 1/14 when the antisymmetry 
of the final two nucleons is taken into account. The ratio becomes 1/5 when the 
antisymmetry is neglected, while the total width is changed only at the level of 20%. 

This discrepancy motivated the search for new mechanisms driving the AA^ — » A^A^ 
decay, such as the inclusion of other mesons, correlated two pion exchange, quark 
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Table 6: OPE results for the ratio r^/Fj; 



TT uncorr. 



vr corr. , 



Takeuchi et al. |5[ (^He) 



Inoue et al. ^ 
Parrefio et al. 
Dubach et al. 



(iHe) 
(iHe) 



(fC) 
(iHe) 

(nuc. matt. 

fl2( 



(iHe) 



0.18 
0.12 
0.09 
0.22 
0.07 



0.06 
0.12 
0.08 
0.10 
0.07 
0.06 
0.20 
0.05 



EXP 



iHe 



iHe 



0.06 ±0.30 
0.25 ±0.13 

0.93 ±0.5 pq 

1 00+1.12 
-'-•'J'J-0.81 
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1.87^ 



-0.91 
-1.59 

0.70 ±0.30 
0.52 ±0.16 
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model calculations, incorporation of Al = 3/2 amplitudes, etc, all of which will be 
described in coming sections. 



3.1.2 Meson exchange model beyond OPE 

The OPE mechanism is an excellent starting point to study the weak AA^ — > A^A^ 
decay due to the experimentally known weak vertex. However, this mechanism can 
at most describe the long range part of the transition potential, while it was soon 



realized ||53| that the large momentum exchanged (400 MeV/c) would probe quite 
short distances and more massive mesons might contribute in the decay process. 

a) The p meson 

Including the p meson in the weak decay mechanism was fist attempted by McKel- 



lar and Gibson |5J|, where the weak vertex Hamiltonian was 



n 



w 

ANp 



' 2M ' ' 



rpf,ipA (?) 



(40) 



where, contrary to the pion, the weak constants a, (3 and e must be determined 
theoretically. Taking the usual strong AA^p vertex 



the nonrelativistic reduction of the p transition potential reads 



(41) 



{a + (3){F, + F^) 

-= ((Ti X g)((T2 X q) 



AMM 



+i 



2M 



[ai X a2)q 



g ^ ± ml 



(42) 
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Using ((?! X g)((T2 x g ) = {aia^) (p — (c?i^)(o^2^)5 and decomposing (c?ig)((T2^) into 
central and tensor operators, this potential contains, as in the case of pion exchange, 
central, tensor and parity violating pieces. 

Two weak coupling models are used in ref. [Q, namely, a modified factorization 
approximation — omitting the factor sin Qc cos Qc — and a pole model which applies 
only to the PC amplitudes. The resulting coupling constants differ substantially in 
both approaches, the pole model giving a combination (a + /5) which differs in sign 
and is a factor of 20 smaller than that obtained in the factorization approximation. 
Combining the OPE results with the p contribution with an arbitrary relative sign, 
gives, for the factorization approach, a non-mesonic rate which ranges between 0.7 
(for 71 — p) and 2.3 (for vr + p). Aware of the limitations of the factorization model, the 
authors of ref. |Q state that their estimates can vary up to a factor 3. An important 
observation was that, similarly as for their OPE results, only the tensor transition 
gave a non-negligible contribution to the p exchange weak amplitude. 

Despite the large uncertainty in the result, an application of this model to finite 
nuclei was performed in ref. where the decay rate of helium hypernuclei was 
explored. The (vr + p) model gave an amplification of their vr-only result and the 
(vr — p) one lead to a complete cancellation. Although the ^He result was in nice 
agreement with the data in the (vr + p) model, the ^He rate was overestimated by a 
factor of 2. 

NarduUi |^ obtains the weak p couplings via a pole model, including, in addition 
to the ground state baryon octet, the first excited 1/2+ baryons. However, it must be 
recalled [B4] that their value for the K ^ tt transition is an order of magnitude larger 
than what is obtained from the kaon decay mode K ^ nn [^]. The PV couplings 
were also calculated within a pole model including the negative parity l/2~ baryon 
poles. By rescaling the results of ref. according to the new coupling constants, 
Nardulli predicts the non-mesonic decay rate of a A in nuclear matter to be 0.7. 

The p meson was also investigated in ref. |57| as a part of their full one meson 
exchange potential, which will be discussed in more detail below. Only the ground 
state baryons are considered in their determination of the PC coupling constants, 
while the PV ones are obtained applying soft meson theorems and current algebra 
techniques. Although the baryon pole diagrams could, in principle, also contribute 
to the PV couplings, they were estimated to be only a several percent of the leading 



current algebra contribution p5| . The contribution of the p meson, acting mainly on 
the S ^ D transition, was found to decrease the pion-only rate by 15%. 

A thorough study of the p- meson contribution to the decay rate of j^C was per- 
formed recently using the model of ref. for the weak ANp couplings and 
employing strong couplings from realistic YN potentials [^,^. One of the findings 
of that work is that the central amplitude of the p-meson contribution cannot be 
neglected and is, in fact, larger than the tensor one. The reason is that the function 
j2{kr), representing the outgoing relative NN wavefuntion, eliminates strength of the 
tensor potential from the short distances. Due to its shorter range, the tensor con- 
tribution of the p is strongly reduced with respect to the central transitions, which 
are governed by the function jo{kr). As can be seen in Table |^ the central rate for 
the p contribution turns out to be 3.5 times the tensor one, once one includes SRC, 
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FF and FSI (considered in a simplified way through the NN correlation function 
/(r) = 1 — joiQcT"), with = 3.93 fm~^). The table shows that the p-meson reduces 



the TT-only rate by 10%, similar to the 15% reduction reported in ref. |£6[, although a 
large fraction of the p-rate came from S ^ D transitions in the latter case. Similar 
qualitative conclusions are obtained when more realistic NN wavef unctions, based on 
a T-matrix calculation employing the Nijmegen93 potential |7B|, are used [Q. There, 
the central to tensor relative strength of the p meson is 1.5 and it reduces the vr-only 
rate by 3%. 



Table 7: The p, tt and tt + p contributions to the decay rate of from ref. []86 |. 





P 

uncorr. 


corr. 


TT 


71 + P 


Central 


0.658 


0.085 


0.004 


0.059 


Tensor 


0.155 


0.024 


0.685 


0.454 


PV 


0.109 


0.011 


0.421 


0.477 


Total 


0.921 


0.119 


1.110 


0.991 



It is clear that other conclusions could be reached had one used a different model 
for the ANp vertex. However, since the central and tensor potentials are weighted 
by the same {a + (3) combination, it is quite likely that both central and tensor 
amplitudes are important and need to be included before any further improvement 
on the potential is made. 

b) The full one meson exchange potential 



N N 



K,K* 



(a) 



(b) 



Figure 9: Non-strange (a) and strange (b) meson exchange contribution to the 
AA^ NN weak transition potential. The weak vertex is indicated by the circle. 



A one meson exchange model to describe the AA^ — > A^A^ transition was first 



reported in ref. |^6|] more than 10 years ago, and the details have recently become 
available As shown in fig. pseudoscalar {71,1], K) and vector {p,uj,K*) meson 
exchanges were included and the A/ = 1/2 rule at the weak vertex was assumed. The 
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PC part of the vertices was obtained from a pole model including both baryon and 
meson pole diagrams shown, respectively, in figs. |10| and |Tl|, where the cross indi- 
cates a weak baryon— >baryon or meson^meson weak transition. Employing SU(6)^ 
symmetry and assuming PCAC, current algebra methods allow one to relate the 
baryon— >^baryon amplitudes to the experimentally known A and E decay p-wave am- 
plitudes, and the meson— i>meson ones to the K ^ n one, which is determined from 
the physical K ^ nn decay rate. The strong couplings for the pseudoscalar mesons 
are obtained from SU(3) symmetry and the Goldberger-Treiman relation, while those 
of the vector mesons require the vector dominance assumption. 



jt,p 



(a) 



jt,p 



(b) 



N N 

11,(0 

N 

)( 

A N 
(c) 



r|,(o 



(d) 



K,K* A 



A N 
(e) 



(f) 



Figure 10: Baryon pole diagrams contributing to the PC weak vertices in the 
AA^ — > A^A^ transition amplitude. 



This model was applied to obtain the non-mesonic decay rate of a A in nuclear 

although 



matter [^,^ and, later, predictions were also given for ^He and 
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no details are yet available. 

Following the lines of refs. 0,|5^ , a full one meson exchange potential has recently 
been constructed to study the non-mesonic decay rate of finite nuclei where the 
nuclear structure details of the hypernucleus under consideration are carefully taken 
into account. A convenient compact expression for the potential is given by 



(43) 



where the index i runs over the different mesons (vr, rj, K, p, u, K*) and a over the spin 
operators Oa € {1 (Central spin independent), aia2 (Central spin dependent), 5*12 
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N 



A 



N 



N 



Figure 11: Meson pole diagram contributing to the AA^ NN transition amplitude. 



(Tensor), c?2^ (PV for pseudoscalar mesons) and (ai x a2)r (PV for vector mesons)}. 
The isospin operator ij^^ is 1 for isoscalar mesons {r],uj), T1T2 for isovector mesons 
(7r,p) and a combination of 1 and T1T2 for isodoublet mesons {K,K*). 

The detailed form of the potential and the explicit values of the coupling constants 
can be found in refs. |^8|,^. To facilitate the comparison, it is worth pointing out 
the differences with respect to the model of ref. [^. In the first place, the strong 
coupling constants employed in ref. are taken from realistic YN interactions, such 

ones. Secondly, the potential of ref. 



as the Nijmegen |87] or the Jiilich [8S 



571 does 



not contain form factors in spite of the large momentum transferred in the process 
which makes the pion to be far off-shell. Thirdly, a stronger correlation function is 
used in ref. which in part can compensate for the lack of form factors, 
work of ref. 



The 



uses a parametrization [^] which lies in between realistic correlation 
functions obtained from G-matrix calculations employing the Nijmegen soft- and 
hard-core potentials. Finally, the contribution of the meson pole diagrams to the PC 
couplings were found to be small and, hence, they were omitted in ref. ||58[| . 

Results obtained from both models are shown in Table |^, were the A decay rate 
in nuclear matter of ref. is compared to that in of ref. |^ for a sequential 
(and arbitrary) addition of the different mesons. Although the reduction with respect 
to the pion-only rate is similar when all the mesons are included (25% in the case of 
ref. [0 and 18% in the case of ref. |0), a stronger sensitivity to the strange mesons 



is found in ref. 58 



The work of ref. includes the mesons by isospin pairs [(vr,p), {K,K*), {ri,uj)] 
and finds important destructive interference effects for each isospin-like pair in the 
tensor transitions, which provide the most important contribution for pseudoscalar 
mesons. On the contrary, the PV rates tend to interfere constructively. The study 
also demonstrates the effect of including realistic A^A^ scattering wave functions in the 
final state, reducing the pion-only rate of Table ^from 1.11 to 0.89. The final results 
of ref. |5^ for the non-mesonic decay rate of are shown in Table |^, using the 
Nijmegen (Jiilich) strong coupling constants. The ordering of the mesons is different 
from the previous table and reflects the isospin-pair study mentioned above. Table U 



33 



Table 8: One meson exchange results for the non-mesonic decay rate. 



Meson 


ref. [El 
(nuc. matt.) 


ref. \m 

aC 


71 


1.85 


1.11 


+7] 


1.68 


1.06 


+K 


1.31 


0.54 


+P 


1.11 


0.51 




1.19 


0.58 


+K* 


1.38 


0.91 



shows the significance of SRC,FF and FSI in the decay rate, drastically reducing the 
influence of the heavier mesons on the rate. The K meson lowers the rate by almost 
50%, reduction that is partly compensated by its isospin partner, the K*. The final 
result is either 15% smaller (using the Nijmegen couplings) or 15% larger (using the 
Jiilich couplings) than the pion-only result. This sensitivity to the strong couplings 
is unfortunate since one expected to use the non-mesonic decay of hypernuclei to 
learn about the weak vertex. More data on YN scattering is highly desirable in order 
to reduce the uncertainties on the strong coupling constants. The table also shows 
the results when new i^-meson weak couplings including one-loop corrections to the 
leading order in Chiral Perturbation Theory are employed. Since the one-loop 
corrections decrease the tree-level p-wave amplitudes by 30%, the influence of the 
/T-meson is reduced substantially giving rise to an increased rate. Note, however, 
that the effect of counterterms was ignored in ref. [^] and further investigations 



are necessary. In any case, any theoretical improvement in the determination of the 
weak and strong vertices can be easily incorporated in the meson exchange model of 



ref. 58 



Table 9: One meson exchange results for the non-mesonic rate of ]^C, obtained in |^ 
using the Nijmegen (Jiilich) strong couplings 



Meson 


uncorr. 


corr. 


TT 


1.68 (1.68) 


0.89 (0.89) 


+P 


2.06 (2.33) 


0.86 (0.83) 


+K 


1.34 (1.70) 


0.50 (0.51) 


+K* 


2.84 (3.82) 


0.76 (0.90) 


+7] 


2.47 (3.82) 


0.68 (0.90) 


+UJ 


2.30 (4.34) 


0.75 (1.02) 


K-couplings 
from xPT 




0.84 (1.10) 



The available results for the non-mesonic weak decay of finite nuclei in the one 
meson exchange model are displayed in Table |T0|. The results for of ref. |]57 



are 



extremely small when compared with their previous nuclear matter results shown in 
Table ^ They are a factor 7 smaller, which seems unreasonable since one may expect 
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reduction factors of at most 1.5 by going from nuclear matter to a finite nucleus 
^^'^ However, the results of the same group quoted in ref. 



79i, shown between 



as C 

brackets in Table [1^, seem to be more consistent with their previous nuclear matter 
results, with the calculations of ref. |Q and with the experimental data. Taking these 
considerations into account we can say that, in general, the non-mesonic rates of finite 
hypernuclei obtained within the one meson exchange model are in fair agreement with 
the measured values. 



Table 10: Non-mesonic decay rates within the meson exchange model of refs. 



and 1^ 





iHe 
TT all mesons 


A ^ 

71 all mesons 


ref. [p^ 
ref. [5S| 


0.9 0.5 
0.5 0.41 


0.5 (2.0 [^]) 0.2 (1.2 [^]) 
0.89 0.75 


EXP: 


0.41 ±0.14 


1.14 ±0.2 
0.89 ±0.18 |7| 



The one meson exchange model has been recently applied to study the non-mesonic 
decay of the hypertriton, [Q. The wave function of the hypertriton, obtained 
by solving the Fadeev equations based on realistic NN forces and the Nijmegen 
soft-core YN interaction, was found to reproduce nicely the A separation energy, 
Ba 



0.13 ± 0.05 MeV I^Tj. The two non-mesonic decay modes, 



d + n and 

nnp, were investigated and the results are shown in Table |Il|. The one-pion only 
result turns out to be 1.7% of the free A decay rate, which is almost a factor 3 smaller 



than what was obtained in a simpler calculation [92], which ignored FSI. Indeed, 
the results in the table show that FSI have a strong infiuence on the rate, reducing 
the PWIA value by a factor of 2. Including all mesons decreases the one-pion result 
by around 10%. That work also makes extensive investigation of the kinematical 
regions where protons and neutrons coming from the partial neutron- and proton- 
induced mechanisms could be found. It was shown that the regions covered by the 
two mechanisms largely overlap and, therefore, the partial rates cannot be measured 
separately unless one applies kinematical detection constraints, which would then 
lead to the measurement of only fractions of these partial rates. 



Table 11: Non-mesonic decay rate of from |^ 





TT 


all mesons 


-pn+d 
^ PWIA 
pn+d 


0.0047 
0.0012 


0.0037 
0.0017 


pn+n+p 
PWIA 
pn+n+p 


0.036 
0.018 


0.028 
0.015 


ptot 


0.019 


0.017 



The construction of a meson exchange model was partly motivated by the unability 
of the OPE mechanism to explain the experimentally large neutron to proton ratio. 
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A summary of results for Tn/Tp is shown in Table 
this ratio, especially for j^^C, obtained in ref. 
is not confirmed by the results of ref, 



57 



The spectacular increase in 
when all the mesons are included 



It could be that the use of form factors, 
considered in ref. |^8| and neglected in ref. |]5^, diminish the effect of the heavier 
mesons, due to the removal of strength at high momentum values and thus at short 
distances, which would then produce a ratio Tn/Tp largely dominated by the OPE 
value. In any case, the difference between the nuclear matter and ^^C results in 
ref. is somewhat surprising in view of the fact that we are dealing with a ratio 
of rates for which the short distances play the most relevant role. Details on the 
finite nucleus calculation of ref. fS^ are awaited in order to clarify the origin of such 
differences. 



Table 12: Neutron- to proton-induced ratio Tn/Tp. 





TT 


7r,p 


all mesons 


nuc. matter 
ref. [57| iHe 

A ^ 


0.06 
0.05 
0.20 


0.08 


0.34 
0.48 
0.83 


ref. 58 1 iHe 

A ^ 


0.104 


0.095 


0.07 
0.07 


EXP: 

m iHe 

12/^ 
A ^ 

[72] i^C 
[B2] i^C 


0.93 ±0.5 

1 qq+1.12 
-'-•'J'J-0.81 
1 Q7+0.91 

0.70 ±0.30 
0.52 ±0.16 



The experimental values for Tn/Tp have large uncertainties due, mainly, to the 
difficulties in detecting neutrons. Less uncertain are the results for the proton-induced 
partial rate Tp and the available data are compared in Table |l^ with the theoretical 
predictions of the meson exchange model of ref. . 



Table 13: Proton-induced ratio Tp. 





EXP 


ref. 5^ 
TT-only all mesons 


iHe 
iHe 
i^B 

12p 
A ^ 


0.16 ± 0.02 g 

0.21 ± 0.07 Q 

0.30l[;:}? |30| 
0.3ll[!-}? 30| 


0.43 
0.64 
0.80 


0.38 
0.56 
0.71 



Obviously, the fact that the non-mesonic rate is well reproduced by this model 
while the ratio Tn/Tp turns out to be too small, leads necessarily to an overestimation 
(of almost a factor of 2) for Tp. It is somewhat surprising that, while both neutron- and 
proton-induced partial rates appear to be in disagreement with the data, their sum 
conspires to give a non-mesonic rate which reproduces the measurements. Although a 
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more proper comparison should be made after the experimental cut offs in the energy 
of the detected particles are implemented in the theoretical models, it is important 
to stress that precise measurements of Tp, more feasible experimentally than the 
ratio Tn/Tp, would bring valuable information to better understand the mechanism 
inducing the AA^ NN transition. 

c) Two pion exchange 

The two pion exchange mechanism of fig. |l2|(a), containing a strong AA^ —>■ SA^ 



transition plus a weak SA^ NN one, was studied by Bando et al. together 
with the usual one pion mechanism for AA^ — A^A^ . The empirical AA^vr and SA^vr 
weak vertices were used and the AI = 1/2 rule enforced. A closure approximation 
was employed to deal with the sum over intermediate states. The potential for the 
AA^ — > EA^ — i> A^A^ transition was obtained for several partial waves and compared 
with the OPE ones. Important interference effecs for the final 1 = channels (only 
accessible to Ap np), constructive for ^Si Si and ^5*1 Pi and destructive 
for ^5*1 Di, were found. As a result the calculated proton-induced rate, Tp, for 
^He is larger than the OPE value in about a factor 1.5. Moreover, the destructive 
interference in the 1 = 1 ^Si — Pi channel produces an overall decrease in the 
neutron-induced rate, r„. Therefore, the A — S coupling reduces the ratio Tn/Tp 
even more, from the OPE value 0.13 to 0.04. 




A 



N N 



N N 




(a) 



(b) 



Figure 12: Two pion exchange amplitudes to the AA^ NN transition through 
coupling to intermediate SA^ (a) or AA^ (b) states. 



In ref. ||5^, this A — S process was considered together with the mechanism AA^ — »• 
A^A —>■ NN of fig. [l^(b). The loop in figs. was computed explicitly and initial 
AA^ and final A^A^ wavef unctions, correlated by the effect of the strong interaction, 
were used. Note that a potential source of doublecounting in the A — S diagram may 
arise because realistic AA^ correlated wave functions have the AA^ EA^ transition 
built inside. The results for the A — E diagram are qualitatively different than those 
of ref. [93| since a reduction of 30% is found in both the ^Si — Si and ^5*1 —^^ Di 
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amplitudes. The A — N ^Sq Sq amplitude is about a factor 1.5 larger than the 
OPE one and would give rise to an increase on r„, although no rates, only amplitudes, 
were calculated in ref. [^ . 

In analogy with recent developments in the NN potential, in which the p and a 
mesons are considered in terms of the exchange of two correlated pions, there have 
been recent attemps to incorporate this picture in the weak AA^ NN transition 
potential The vrvr correlation is approximated, as displayed in fig. |13|, by a 



t-channel resonance which can be either the scalar a meson or the vector p meson. 
The weak vertex is obtained in both cases by means of a loop diagram with two pions 
and an intermediate or S baryon. The strong vertex is generated in the same 
way in ref. [Q, while phenomenological crNN and pNN couplings from ordinary one 
meson exchange potentials that fit the A^A^ scattering data are considered in ref. [|6T |. 
Note that the strength of the 27r contribution obviously depends on the cut off used 
to regularize the integrals in the triangle loops. 



N 



N 



Figure 13: Correlated two pion exchange AA^ 
couple to 0" or p mesons. 



A^A^ transition, where the pions 



In ref. 60 



a substantial contribution of the 27i diagram to the ^5*0 — >^ 5*0 ampli- 
tude, about 4 times larger than the OPE one, is found. For the other channels, the 
27r contribution represents at most 30% of the vr-only value. This is in contrast to 
what is found in ref. where the 27r/ p potential in the ^5*1 Di channel tends to 



cancel quite strongly the vr-only contribution. The ^5*0 — Sq transition is also quite 
affected by the 27i/p and 27r/(T contributions, which have opposite sign than that for 
the TT-only one, but not at the level of the factor 4 found in ref. ^0 . 
The model of ref. 



ofi^ 



|61 



is applied to compute the partial non-mesonic decay rates 
C It is found that the 27r/p amplitudes decrease the non-mesonic rate from 
0.86 to 0.79 and including also the 27r/cr ones gives a final value of 0.92. The ratio 
Tn/Tp goes from the value 0.087 (vr only) to 0.10 (tt + 2n/p) and to a final 0.14 
(tt + 27r/p + 27r/cr), still far away from the empirical results. 
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3.1.3 Quark model based results 



In recent works |^^^] an old idea of treating the AA^ 
point of view of quark degrees of freedom 



63 



A^A^ mechanism from the 
is reexamined. In the latter reference 



a hybrid quark-hadron model was used in which the usual OPE between hadrons 
was adopted for the external region of relative AA^ distances (r > rg with tq — 0.8 
fm), while a six-quark wave function and an effective quark Hamiltonian is taken 
for To < 0.8 fm. The main uncertainty within that model is the parametrization of 
the weak effective AS = 1 quark Hamiltonian. Only when the coefficients of the 
QCD corrected Oilman and Wise Hamiltonian are modified to implement the 
A/ = 1/2 rule, did the non-mesonic decay of fall within the uncertainties of the 
experimental Brookhaven value of 1.14 ± 0.2 . 

In refs. p4| , |66[| a four-quark AS* = 1 effective Hamiltonian, containing both A/ = 
1/2 and 3/2 terms and modified by QCD effects [^] was used, along with the quark 
cluster model for the quark component of the two-baryon initial and final states. Since 
quark- ant iquark pairs are not allowed as intermediate states, the authors incorporate 
explicitly meson exchange diagrams, as the vr meson and also the K meson . 

Both works find large A/ = 3/2 contributions to the ^5*0 — >^ 5*0 and ^Sq Pq 
amplitudes. However, it would be desirable to clarify the origin of some sign and size 
differences for certain amplitudes. For instance, the neutron-induced contribution to 
the ^5*0 -^^ 5*0 amplitude is —0.07 times the proton-induced one in ref. |^^, while 
it is a factor —0.75 the proton-induced one in ref. |^. The proton- induced ^S'o — > 



Pq amplitude is much larger in size in ref. \EB\ than the one obtained in ref. [p^|. 
Therefore, although the results of both works are qualitatively consistent, they still 
show important quantitative differences. 



Table 14: Quark-model results of ref. |£5|] for the non-mesonic decay rate and ratio 
Tn/Tp of light hypernuclei. 





OPE 


OPE and DQ 


EXP 


rn?Tt(AHe) 

rn/r,(iHe) 


0.216 
0.132 


0.627 
0.489 


0.41 ±0.14 
0.93 ±0.55 


CO c3 


rnm(AHe) 

rn/r,(iHe) 


0.154 
0.061 


0.253 
0.178 


0.20 ±0.05 
0.25 ±0.13 


81] 0.17 ± 0.05 g 

81] 0.06 ± 0.30 [8C] 



The amplitudes obtained in ref. [34 



were used to determine the partial non- 
mesonic decay widths of s-shell nuclei but the results depended on the relative sign 
between the quark and the OPE contributions. Using the soft pion theorem, the 
authors clarified recently ^5[ the sign ambiguity. This theorem allows one to relate 
the sign of the A 



rnr 



coupling constant (relevant for the OPE mechanism) to 
the PC baryon-baryon transition amplitude A — >■ n, which was computed using the 
quark Hamiltonian and their quark model wave functions. A more realistic A wave 
function, pushed to the surface by the effect of a repulsive A-nucleus potential at short 
distances, was also used in ref. [Q. The results, summarized in Table 0, improve 
substantially those obtained in their previous calculations These results are 

promising but, before calculations of the decay of heavier hypernuclei are attempted 
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within this model, it would be desirable to establish the connection between the 
effective quark Hamiltonian and the empirical A Ntt vertex. 



3.1.4 A/ = 1/2 violation? 



a) Phenomenology 

Non-mesonic weak decays of light hypernuclei can be used to test the validity of 
the phenomenological A/ = 1/2 rule. One must first observe that, since the AA^ pair 
is in a L = state, the possible transitions are those listed in Table [15. 



Table 15: Possible AA^ NN transitions starting from a L = AA^ pair. 



'So - 


^Sq 


(1=1) 






(1=1) 


'S, - 


- 'S, 


(1=0) 




- 'D, 


(1=0) 




- 'Pi 


(1=0) 




- 'Pi 


(1=1) 



The neutron-induced transitions can only lead to / = 1 final states and, if the 
A/ = 1/2 rule applies, then the ratio of neutron- to proton-stimulated transitions to 



final 1=1 states must be 52 



p(7=l) 

7^ = ^- (44) 
i p 

Therefore, by isolating the initial singlet spin state ^5*0, which leads to only 1 = 1 
final states, both in neutron- and proton- induced decays, the validity of the A J = 1/2 
rule in the non-mesonic decay could be verified. This idea stimulated a large amount of 
experimental effort pO| , |8T| , p6| , [97[| . The connection with experiments is done following 
the phenomenological model of Block and Dalitz [52| and updated by Dover ^8 



according to which the non-mesonic decay rates of light hypernuclei are expressed in 
terms of elementary interaction strengths between singlet or triplet AA^ pairs: 

r„™(iH) = lp,{3Rni+ RnO +2Rpo) (45) 

o 

r,^(iHe) = Ip3{ 2Rno + 3Rpi+ Rpo) (46) 

o 

rnm(AHe) = -p4,{3Rnl + RnO + 3-Rpi + Rpo) ■ (47) 
o 

The quantities Rns stand for the non-mesonic decay of a AA^ pair with spin S, for unit 
density of nucleon A^ at the A position, and pA denotes the mean nucleon density at 
the A position. An average over spin and charge is implied in eqs. (|45|) — d^) . These 
expressions also assume that the final state interactions, which were shown to be 
important in the non-mesonic decay of the hypertriton |9^, are contained in the 
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elementary interaction strengths and thus are the same in all these light hj^pernuclei. 
The densities cancel in forming the ratios 



y ^ r„(iHe) _ 2R^o ^^g^ 
Tp{\}ie) SRpi + Rpo 

^5 ^ r„(|He) ^ 3Rnl + -RnO ^^gsj 

Fp^He) 3i?pi + -Rpo 
^ _ TrimdHe) _ 2i?„o + 3-Rpi + Rpo ^^^^ 

^nmiA^) 3i?„i + RnO + 2i?p0 



all of which involve experimentally known decay rates, as seen in Table |16|, and 
therefore allow one to extract the ratio 

RnO Inml^ /c-tn 

K = = . (51) 

RpO 1 + 7^ - 7nm7 



Table 16: Experimental weak decay observables of several s-shell hypernuclei. 





iHe 

m 


4 

H ^ m 


iH 

i3 


r 

nm 


0.20 ±0.11 
0.21 ±0.07 
0.93 ±0.55 
0.41 ±0.14 


0.04 ±0.02 0.01 ±0.05 
0.16 ±0.02 0.16 ±0.02 
0.25 ±0.13 0.06 ±0.30 
0.20 ±0.05 0.17 ±0.05 


0.17±0.11 



Taking the central values of the KEK ||80| or BNL experiments, the value for 
K is 0.4 or 0.8 but the error bars are so large that a value of 2, compatible with the 
AI =1/2 rule, is not ruled out |97]. A much more direct determination of k would 
be obtained from 

I He) RnO 



rp(iH) 



R 



'pO 



(52) 



which requires the measurement of the partial rates for ^H, an attempt that will be 



made in experiment E907 at BNL [^. Since the partial neutron rate for ^He has 
been measured to be small, even if the error bar is large, the A J = 1/2 rule would 
pose severe constraints on the value of Tpijji), which could not exceed a few percent 
of the free A decay rate. 

b) Models 

The quark model of refs. |^,^ is based on the effective 4 quark Hamiltonian 
with QCD corrections |^| extended down to an energy scale /i^ ~ /Iq |95] at which 
«(/iQ) = 1. This quark Hamiltonian contains some enhancement of the AI =1/2 
transitions although it is still insufficient to explain the ratio of AI = 1/2 to AI = 3/2 
amplitudes for the non-leptonic weak decay of K, A and other strange hadrons. As 
mentioned in the previous section, two body amplitudes were evaluated by combining 
the quark Hamiltonian with a constituent quark model for the baryons. In spite of 
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the numerical differences, both models predict substantial A/ = 3/2 contributions to 
the spin singlet amplitudes ^5*0 5*0 and ^5*0 — Pq. Therefore, these are the most 
important ones to measure in order to determine violations of the A/ = 1/2 rule. 
In ref. |]9^ violation of the A/ = 1/2 rule in the AN — > NN mechanism was 



investigated in the framework of a meson exchange model. The work is based on the 
observation that, while the A/ = 3/2 contributions to the weak coupling constants 
of baryons to the pseudoscalar mesons are empirically known to be small, this is 
not necessarily the case for the vector mesons. Indeed, based on the weak QCD- 
corrected four quark effective Hamiltonian at a scale ~ 1 GeV [|100|| , the factorization 
contributions to the weak ANp and TiNp couplings were found to contain A/ = 
3/2 terms comparable in magnitude to the A/ = 1/2 ones [Q. The factorization 
approximation in the evaluation of the A/ = 3/2 amplitudes for the known hyperon 
decays yields |p.01|| : (1) good fits for the s-wave and p-wave S amplitudes and the s- 
wave triangle discrepancy, (2) an underestimate of the p-wave E triangle discrepancy 
by a factor 3-4, and (3) an overestimate of the s-wave and p-wave A amplitudes by 
a factor 3-4. It is important to notice, however, that the experimental errors on the 
A/ = 3/2 amplitudes are large and the above mentioned discrepancies are within 



the 2a level. The work of ref. |99| takes the predicted A/ = 3/2 contributions to 
the coupling constants of the vector mesons and rescales them by factors ranging 
in between (—3, 3) to account for the limitations of the factorization model. Their 
A/ = 1/2 potential was extended to incorporate the A/ = 3/2 transitions of the 
p and K* mesons by assuming the A to behave asa|3/2 — l/2> isospurion and 
replacing the isospin rr operator in the potential by 73/2^^, where 73/2 is the usual 
isospin 3/2 — > 1/2 transition operator. Even if the A J = 3/2 couplings of the 
vector mesons were comparable in magnitude to the corresponding A J = 1/2 ones, it 
was found that, when all the other mesons were incorporated, the decay rate of 
changed at most by 6%, the proton-induced rate was affected at most by 10%, while 
the neutron-induced rate could change by up to a factor 2. In the most favorable 
scaling situation (a scale factor 3 for both p and K*) the ratio F^/Fp changed from 
the A/ = 1/2 value of 0.07 to 0.14, which is still away from the experimental values. 

These findings are substantially different than those found in the quark model 
works 0,0], were the A/ = 3/2 amplitudes seem to play a much more relevant 
role. Therefore a direct determination of the spin singlet ratio Rno/Rpo, through the 
measurement of nucleon induced decays of light hypernuclei, is very much desired in 
order to clarify the issue of A/ = 1/2 violation in the non-mesonic weak decay. 



3.2 Asymmetry 

The possibility of measuring the asymmetry in the angular distribution of protons 
from the weak decay of polarized hypernuclei offers a unique opportunity to learn 
about the weak decay mechanism. The asymmetry is due to an interference between 
PC and PV amplitudes and, hence, can give complementary information to that 
obtained from total and partial rates which are dominated by the PC piece of the 
AA^ NN interaction. 

The first observation of asymmetries on the decay products from polarized hyper- 



42 



nuclei was done at KEK |7^. There, polarized hypernuclei (^^C, ^ B) were created 
through the (tt"*", K~^) reaction, which has been shown to produce nonnegligible po- 
larization perpendicular to the reaction plane for a momentum of the pion p^^ ~ 1.04 
GeV/c at small K'^ scattering angles (10 — 15°) where the cross section is still appre- 
ciable |T02| , M . 

The intensity of protons emitted at an angle x with respect to the polarization 
axis in the decay of polarized hypernuclei is given by 

lix) = TiiMpM^) , (53) 
where Ai is the hypernuclear transition operator and 



is the density matrix describing the hypernucleus of spin J polarized along the y axis, 
which is perpendicular to the reaction plane. Eq. (^) can be written as 

Iix) = Ioil + PyAyix)) , (55) 

where Iq is the isotropic intensity of protons from the decay of unpolarized hypernuclei 
and Ay{x) the asymmetry parameter which reads 

y a(Mi)Mi 

^^^^^"j+1 MMM^) "J+1 ya{M,) ^56) 

The last equation shows that the asymmetry exhibits a simple cosx dependence 
and is given in terms of a{Mi), the weak decay probability for a hypernucleus with 
spin projection Mi ejecting protons along the quantization axis. The experiment 
compares the total number of protons emerging parallel or antiparallel to the y-axis 
and , therefore, the product PyAy is determined. Using standard angular momentum 
algebra and assuming the A to couple only to the nuclear core ground state of spin 
Jc, one can obtain the A polarization as 

PA = { ~jTT^' ^ = " (57) 

Py if J= Jc + 1/2 

and define the intrinsic A asymmetry parameter 

a,= l-'^^y ifJ=Je-l/2 (58) 
[ Ay if J = Jc + 1/2 



such that the asymmetry 



A = PyAy cos X = Pa^a cos X (59) 
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is expressed in terms of A properties, pa and oa, which are then characteristic of the 
elementary AA^ — > A^A^ reaction. 

In order to compare with experiment, the polarization of the hypernucleus must 
be known. Until now this information is not available and it is necessary to resort to 
theoretical models, which take into account the conditions of the production reaction. 
Given the energy resolution of the KEK experiment (5 — 7 MeV), several hypernu- 
clear states can be excited, which then decay electromagnetically and/or by particle 
emission to a hypernuclear ground state prior to the weak decay. To determine the 
polarization at this stage one requires: i) the amount of polarization of the states ex- 
cited by the production mechanism, along with their corresponding formation cross 
sections, and ii) an attenuation coefficient to account for the loss of polarization of 
each level in the deexcitation process [|104|] . 

In ref. |102| one-particle one- hole (j~^Ja) hypernuclear wave functions were em- 



ployed to obtain the polarizations of typical states within the p-shell region. The 
model was refined in [|103|| by incorporating configuration-mixed shell-model wave 
functions. Using their predictions for the three 1~ states in and the depolar- 
ization formalism of ref. ||104| , one obtains Py = —0.19 {pa = 0.095). On the other 
hand, the hypernucleus is created by particle emission from excited states of ]^C. 
The window of excitation energy that spans 1.55 MeV between the (p +a^B) and the 
(A+^^C) particle decay threshold contains three positive-parity states: two 2+ states 
separated by ~ 800 KeV and a narrow 0^ state just below the (A+^^C) threshold. 
Using the model of ref. |[L03| ], which predicts equal formation cross sections for the 
2f and 22 states, and neglecting the O"*" state because of its relatively small cross 
section, a polarization of Py = —0.29 {p\ = 0.21) is obtained. However, hypernuclear 
structure calculations by Auerbach et al. ||105|| predicted strong configuration mixing 
which reduced the cross section of the lower state by a factor of three relative 
to the higher one. This prediction was verified by a reanalysis of older emulsion 
data [|106|| . Taking these relative weights into account, one obtains the more realistic 
value Py = -0.43 {pa = 0.31). 

In Table |l^ a comparison of theoretical results to the presently available data is 
made. The two meson exchange models predict a quite similar asymmetry parameter 
a\ which ranges between —0.3 and —0.4 for the different hypernuclei. The comparison 
with experiment must be made at the level of the asymmetry, for which it is necessary 
to multiply by the A polarization discussed above. Good agreement with the data 
is obtained. However, the error bars are still so large that not much can be learned 
about the mechanism governing the decay. 

In order to avoid the need for theoretical input and access Ay directly, a new ex- 
periment at KEK ||107|| was devised to measure the decay of polarized ^He, extracting 
both the pion asymmetry from the mesonic channel, At^- , and the proton asymmetry 
from the non-mesonic decay, A. The asymmetry parameter a^- of the pionic channel 
has been estimated to be very similar to that of the free A decay [|108|| , and, there- 
fore, the hypernuclear polarization can be obtained from the relation Py = A-j^-fa^-. 
Combining this value with the measured proton asymmetry A, permits to determine 
the asymmetry parameter for the non-mesonic decay from the equality Ay = A/Py. 

A preliminar analysis of this experiment has recently become available ||109|| and 
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Table 17: Proton asymmetry from the weak decay of hypernuclei. 





Parreno et al. [B8|] 
iHe i^B 


Dubach et al. ^| 
nuc. matter 


ttA (tt corr.) 
oa (all mesons) 


-0.238 

-0.273 -0.391 -0.316 


-0.192 
-0.443 


Pk 


0.31 0.095 




A = PAttA 


-0.12 -0.13 




A (EXP) |7| 


-0.20 ±0.10 -0.01 ±0.10 





gives: = 0.217 ± 0.087 ± 0.021 for an emerging angle 2 <| 9k \< 7° and 
= 0.382 ± 0.114 ± 0.013 for 7 <| 6k \< 15°. For the proton asymmetry the 
measured values are: A = 0.102 ± 0.050 ± 0.005 for 2 <| |< 7° and A = 0.045 ± 
0.099 ±0.021 for 7 <| 6k \< 15°. The combination of these preliminary results, which 
show large uncertainties, would give rise to a positive asymmetry parameter a\, in 
total contradiction with their previous results |]73[ and the theoretical predictions 
shown in Table |l^. One would also expect a straight line, passing through the origin, 
when representing the asymmetry A = Pa^a as a function of p\. The two new 
experimental points do not even show this behavior. However, one should keep in 
mind that corrections due to the leaking of pions into the proton identification area 
still need to be made 



lOS] and, therefore, the analysis is not yet complete. 



3.3 Two-nucleon induced decay 

The relevance of the two-nucleon induced channel, ANN NNN, in the decay of 
hypernuclei was first pointed out by Alberico et al. [|110|| . This transition can be 



viewed as coming from the absorption of a virtual pion emitted in the weak vertex by 



two strongly correlated nucleons, as depicted schematically in fig. |1J. An interesting 
point raised in ref. ||110[| was that this channel affects the emission of neutrons and 
could have some relevance in the experimental determination of the ratio r^/Fp. Since 
pions are absorbed mainly by neutron-proton pairs, the 2A^-induced decay proceeds 
basically through the mode Anp — > nnp, emitting twice as many neutrons than 
protons. Hence, an experimental observation of a large number of neutrons could still 
be compatible with a smaller value for F„/Fp if the 2A^ mechanism was appreciably 
large. 



In ref. ||110|| the decay width of a A in nuclear matter was calculated within the 
propagator method formalism discussed in Sect. 2.1.1, which allows treating all the 
decay channels in an unified way. Eq. (10) needs to be modified to include the effect 
of form factors and short-range correlations, which have been shown to be important 
for the lA^-induced mechanism AA^ NN . Thus, a monopole form factor at both 
the strong and weak vertices 

^(^) = ' (60) 

with A ~ 1.2 GeV is included. Moreover, the pion interaction between two nucleons 
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A 




Figure 14: Schematic representation of the A decay coupHng to 2p2h components 
through virtual (close to real) pion absorption. 



is replaced by the particle- hole interaction of refs. [|lllj , |112( 

Vph = {VL{q)qiqj + VT{q){Sij - qiqj)} aiajTT , 



(61) 



which includes tt and p exchange modulated by the effect of short range correlations. 
Detailed expressions of Viiq) and Vriq) can be found in ref. [T^. It is convenient 
to split VL{q) into the free pion propagator plus the Landau-Migdal parameter, g', 
which here is slightly momentum dependent 



VL{q) 



fF\q)D,{q)+g\q) 



(62) 



The AA^ correlations are treated in a similar way as in |jT3| by modifying the pion 
exchanged between the weak and strong vertices. The decay width can be finally 
written as 



m 



2\2 



<fq 
(2^ 



[1 - n{k - q)]e{U^ - E{k -q)-Vi 



X Ima(g)|^o=fcO_s(fc_g.)_v^ ' 



(63) 



with 



a{q) = \S'+ -] q'' \ F\q)D,{q) 



+ 



7^. 
~S\q)U{q) 



+ 



Pl{q)U{q) 



+ 2 



1 - VL{q)U{q) 1 - VL{q)U{q) 
PHqWiq) 



1 - VT{q)U{q) 



(64) 
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where the functions S{q), Pl{.(i) and Priq), defined exphcitly in ref. [|T3l, reduce to 



S{q) ^ (^i^yDoiq)F\q) (65) 

pLiq) - (-]-fDoiq)F'iq) (66) 
Pxiq) ^ (67) 

in the absence of AA^ short range correlations. 

The 2A^-induced A decay mode appears when one considers the absorption of the 
virtual pion by 2 particle-2 hole {2p2h) states. Therefore, the function U{q) in eq. ( |64|) 
must include a term, U2p2h{q), that accounts for this coupling and becomes 

Uiq) = U^iq) + U^iq) + U^Mq) ■ (68) 

Before going into the details on how this function is determined, let us illustrate 
the physical mechanism that leads to the new decay channel. The strength of a free 
pion, which is accumulated in a delta function around a momentum value of 100 
MeV/c, spreads in the medium over a wide range of energies due to the coupling of 
the pion to ph, Ah and 2p2h excitations. This distribution presents a Breit-Wigner 
type peak, the position of which changes slightly from the original pion pole. The 
energy and momentum of the pion at the peak is such that, at normal nuclear matter 
density, still enforces a nucleon momentum which lies below the Fermi momentum. 
However, the width of the distribution is such that part of its tail corresponds to a 
Pauli unblocked situation and, since at low pion energies this width is mostly due to 



pion absorption through 2p2h states, as shown in fig. |14|, the new mode would be 
observed as three particle emission from ANN —>■ NNN. 



In ref. ||110|| the following parametrization was used 

U2p2h = -47rg V'Co , (69) 
with a value ImCo = 0.18/i~^ taken from theoretical studies performed in connection 



with (e, e') reactions ||113|| and compatible with data from pionic atoms, provided the 
real part of Cq was ReCo = OAOfi~^. The value of U2p2h was taken constant, even if 
the phase space in the A decay forces one to move away from the situation in pionic 
atoms to values of {q^,q) different from (/i, 0). Moreover, lmU2p2h was set to zero 
in the region where Imf/^r 7^ 0. Their results for the 2A^-induced decay rate of a A 
in nuclear matter at an average nuclear density p = 0.69po; where po is the normal 
density, are shown in Table |1^, splitted into the s-wave, p-wave longitudinal and p- 
wave transverse contributions. The final result of 0.52 represents a fraction of about 
40% of the lA^-induced decay rate. 

The relative large strength of this channel and the repercussion it has in the proper 



interpretation of the ratio Tn/Tp motivated a more detailed analysis [^, in which 
a more realistic pion self-energy was used and the results were extended to finite 
nuclei via the LDA. The determination of U2p2h was based on new developments in 
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Table 18: 2A^-induced contribution to the A decay rate in nuclear mattter. 





m 




Ph. space 


s 


0.40 


0.119 


0.148 


Pl 


0.079 


0.023 


0.026 


Pt 


0.040 


0.023 


0.065 


^2p2h 


0.52 


0.165 


0.238 



the extraction of the optical potential from pionic atoms ||114|| , reanalyzed in ref . ||115 



to account for different neutron and proton radii, which yielded 

1 



n2p2h(g° = /X, g ~ 0, p) = -47rg V' (l + Co (70) 



with 



M 

Co = (0.073 + i0.068)/i-^ 

However, since the propagator method model described above generates the Lorentz- 
Lorenz correction automatically, one needs to rewrite the information of eq. ( ffOD 
accordingly. Thus, for the average density of pionic atoms in p-wave (0.75po according 
to | 116|| ), the results of our formulation should coincide with the experimental one 



given in eq. (^). This is readily obtained with the help of eq. (8) and the result for 
U2p2h is i3 ^ 

g"f/2p2h(g° = /i, g ~ O, p) = -47rg V'C* , (71) 

with 

C* = (0.105 + i0.096)p~^ , (72) 

which yields a value for Im Cq half of that of ref. ||1 10|] and a value for Re Cq four times 
smaller. The extension of the self-energy to new kinematical regions away from pionic 
atoms was also performed in ref. [0. The importance of all these modifications is 



displayed in Table |I8|, where the third column corresponds to the same calculation 
of ref. ||110|| but replacing their Co value by the new Cq one. The results are reduced 
by a factor 3 due to a smaller absorptive ImQ and a less attractive Ren2p2?t which 
moves the pion pole to a smaller momentum. Not forcing lmU2p2h = in the ph 
excitation region (where ImUN 7^ 0) and taking the phase-space modifications into 



account increase the results further, as seen in the last column of Table 0, but the 
final value is less than half that obtained previously ||110|] . 

Even if this decay channel turns out to be considerably smaller than claimed 
before, it will be shown in the next section that its consideration has repercussions 
in the experimental determination of the F^/Fp ratio. This investigation requires 
to have predictions for the 2A^-induced decay rate of finite nuclei, which were also 



obtained in ref. |67] on the basis of the LDA shown in eq. (13). The results are 
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Table 19: Mesonic, lA^-induced and 2A^-induced decay rates for several hypernuclei 
in the IDA 





F 

771 


r 


J- O^nO/i 




0.31 


1.45 


0.27 


i6n 

A ^ 


0.24 


1.54 


0.29 


i^Ne 


0.14 


1.60 


0.32 




0.03 


1.76 


0.32 


i^Fe 


0.01 


1.82 


0.32 


89v 
A ^ 




1.88 


0.31 


208pb 




1.93 


0.30 



displayed in Table W for different hypernuclei, where, for completeness, the decay 
rates from the mesonic and liV-induced channels are also shown. 

The total lifetimes — inverse of the total decay widths — are somewhat small 
compared with the recent lifetime measurements at KEK ||117|| , which in the case 



of ^ Fe is 215 ± 14 ps, a factor 1.8 larger than the result obtained from Table 



However, the lifetime of ^°'^Pb derived from the table is only 20% smaller than the 
central value 145 ± 7 ± 23 ps obtained for ^°^Bi in Jiilich ||118|| and compatible within 
errors. The latter measurements are done by bombarding heavy nuclei with protons 
and looking at the delayed fission fragments. The method was proved to work in 
ref. [ [119|| and used in the experiment of ref. ||118|| in order to obtain the hypernuclear 
widths with more precission. These latter experimental results are also compatible 
with those of ref. [|120|| obtained from p interactions with nuclei which, however, 
have large uncertainties. One must note that the theoretical results are somewhat 
dependent on the AA^ correlations, controlled by the g'j^ parameter, and on the A wave 
function used in the LDA. The lA^-induced rate for would be reduced from 1.45 
to 1.26 if a more extended A wave function with an oscillator parameter 6a = 1.87 
fm and a value g'j^ = 0.2 were used. This value for g'j^ is what corresponds to a Bessel 
type correlation function /(r) = 1 — joiqc"^) with Qc = 3.93 fm~^, very similar to the 
correlation function used in ref. [|58| based on realistic AA^ wavefunctions obtained 



from the Nijmegen interactions. Therefore, the results in Table |T9| could easily be 
10-15% smaller and, correspondingly, the lifetimes 10-15% larger, which would then 
agree better with the very heavy hypernuclei results ||118| , |120|] but would still be in 
disagreement with the KEK ones. 



3.4 The Tn/Tp puzzle 

As mentioned in the previous section, one of the interesting consequences of the 2A^- 
induced mechanism is that, since pions are absorbed predominantly by np pairs, this 
decay mechanism would emit two neutrons and one proton. Hence, part of the large 
number of observed neutrons would then come from the 2A^-induced channel and the 
ratio Tn/Tp would not need to be so large as what is found in the present analyses, 
which associate all non-mesonic decay to the lA^-induced mechanisms. Actually, as 



shown in ref. |^ and summarized below, the situation is not so clear. 
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The experimental analysis of ref. |^ takes only the lA^- mechanism into account 
and the assignment of the number of neutrons and protons goes as 

iVp ~ r;^p (73) 

with = r^^P + Fp^P. However, considering the 2 A^- induced mechanism and as- 
suming that all emitted nucleons are detected, the appropriate analysis reads 

Nn ~ 2r„ + Fp + 2T2p2h 

Nj, ~ Tj, + T2p2h, (74) 

where the 2A^-induced decay rate T2p2h is assumed to proceed predominantly through 
the Anp nnp channel. 

On combining the sets of eqs. (^) and (|7^ one finds 




n \ \ \ 



2 1 1 pexp I I pGxp 



2p2h 



2p2h 

(75) 



pexp I I pcxp 

which gives the new ratio in terms of the old one and the 2iV-induced decay rate. From 

pexp 

eq. ( |75D one can see that if — ^ > 0.5, as is the case of almost all experimental results, 

Fp 

then the ratio F„/Fp of the new analysis is even larger, in stronger disagreement with 
the present theories. Taking as an example the data of ref. ^0[ for ]^C, corrected 

pexp 

for final state interactions ||121|] , one extracts — ^ = 1.04 while the new analysis 

Fp 

gives — = 1.54 which shows that, even if these values have big error bars tied to 
Fp 

the uncertainties in and Np, the 2A^-induced decay channel is relevant and needs 
to be considered in any experimental analysis from where the ratio F„/Fp is to be 
extracted. 

Note, however, that not all the three nucleons emitted in the 2A^-induced decay 
are necessarily fast ones | ]122| |. If the process receives most contributions when the 
pion is emitted close to on shell, then the nucleon at the A vertex will be slow and 
will not be detected due to the experimental energy tresholds of around 30 — 40 MeV. 
Assuming this to be the case, the experimental reanalysis gives, instead of the relation 
in eq. (|75D , the following one 





(76) 



p pexp pexp 

with Tnm = F„ + Fp + F2p2h, which yields < -i^p as long as < 2. 

i p i p i p 
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These considerations lead to the conclusion that, in order to establish the effect 
of the 2A^-induced decay mode on the ratio r^/Fp, sensitive to the detection thresh- 
olds, it is necessary to know the spectrum of the emitted nucleons from the different 
mechanisms. This is discussed in the next section. 



3.5 Nucleon spectra 

The findings of the preceding section evidence the need for theoretical calculations of 
the energy distribution of nucleons from the different weak decay mechanisms. These 
models will also have to address the problem of final state interactions of the nucleons 
on their way out of the nucleus. 

The spectra of neutrons and protons from the weak decay of hypernuclei have 



been calculated in a recent work ||123|| , taking F^/Fp as a free parameter in order to 
facilitate the extraction of this ratio from the experimental data. The spectra were 
obtained using the propagator method formalism |1^,|6^ to determine the primary 
emitted nucleons for a given hypernucleus via the LDA, an approximation particularly 
suited to treat the final state interactions, as we shall see. Let us recall that the rate 
for a particular channel (z = n, p, 2N) is given by 

F, = / d'kp{k)Uk) = I d'kp{k) I d\ I 0A(r) P F,(fc,p(r)) , (77) 



where Fj(A;, p(r)) is the decay rate of a A of momentum k in nuclear matter of density 
p(r), which is the nuclear density at point r where the decay occurs. The function 
0a(^ is the A wave function at that point and an average over the A momentum 
distribution is also implied in eq. (|77D. The structure of the integrals involved is then 

Ti = j d^k J d\ Jd^q...lm Ui{q^ = k^ - E{k - q) - V^v, q) (78) 

where the function Ui associated to the particular decay mechanism can be further 
decomposed in terms of the momentum variables of the emitted nucleons, as explained 
in ref. ||123|| . Then, a Monte Carlo integration method can be used to perform the 



integrals, where the random configurations are generated such that: i) give the relative 
amount of one-nucleon- and two-nucleon-induced decay probabilities according to the 
model described in Sect. 3.3, and ii) give the relative amount of neutron- and proton- 
induced decay rates, according to the input value of F„/Fp, used as a parameter of 
the calculation. 

In this way, each configuration determines, at each point r, a number of primary 
emitted nucleons and their corresponding momenta. The fate of these nucleons on 
their way out of the nucleus is then followed using a Monte Carlo simulation method 



77| , which allows the nucleons to undergo collisions with other nucleons according to 
NN cross sections modified by Pauli blocking. In each collision, the nucleons change 
energy, direction and, eventually, charge, as well as exciting other nucleons above the 
local Fermi momentum. Then, the propagation of these secondary unbound nucleons 
through the nucleus is also followed. 
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Figure 15: Spectra of neutrons and protons in the decay of j^^C, assuming a ratio 
Tn/rp = 1- Dashed hne: protons from lA^-induced decay. Dotted hne: neutrons 
from lA^-induced decay. Sohd hne: protons from 2A^-induced decay. Dash-dotted 
hne: neutrons from 2A^-induced decay. 



The obtained spectra of neutrons and protons coming from the lA^- and 2N- 
induced mechanisms in the decay of are shown in fig. |l^ for a value F^/Fp = 1. 
One observes that the distribution of protons (and neutrons) from the liV-induced 
mechanism peaks around 70 — 75 MeV, which corresponds to the energy of the most 
probable kinematics in which the two nucleons emerge back-to-back. The energy 
distribution generated by to the 2iV-induced process is flatter and presents a peak 
at low excitation energies, which contains mainly slow nucleons coming from the A 
vertex plus some nucleons whose energy has been degraded by the effect of final state 
interactions. It is clear then that the experimental thresholds (30 — 40 MeV) will 
prevent these 2A^-induced nucleons from being detected. However, a considerable 
fraction still leaks into the region where most liV-induced decay nucleons appear. 
Hence, separating the nucleons from both mechanisms will be very difficult, unless 
angular correlation measurements are also conducted. This effort is presently being 
carried out at Brookhaven and it is certainly encouraging to see that their preliminary 
results [81| already show how the back-to-back kinematics nicely selects the nucleons 
emitted in the lA^-mechanism and, therefore, a cleaner extraction of the ratio F„/Fp 
might soon become available. 



Qualitatively similar spectrum shapes were obtained for ^He in ref. [p.24|| , where 



the 2A^-induced mechanism was assumed to proceed via a A-excited intermediate 
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state, although the medium modification of the A and other important 2p2h mecha- 
nisms were ignored. Final state interactions were not considered, which is reasonable 
for the light hypernucleus treated but does not allow one to apply the model to heavier 
hypernuclei. 

It is clear that the spectra shown in fig. [l^ would certainly differ if other values 
of the ratio r^/Fp were used. A comparison of the available experimental proton 
spectra p0|,|8^ with the calculated ones for several values of F„/Fp was also done 
in ref. | |12ij| ] and it was found that values of F„/Fp ~ 2 — 3 were favored, in strong 
contradiction with the OPE models. 

However, in the same work, other observables that were less sensitive to the specific 
details of the Monte Carlo simulation determining the final shape of the spectrum 
were also studied as functions of F„/Fp, which would probably lead to a more reliable 
determination of this ratio. By integrating over the energy spectra, the total number 
of neutrons, A^„, and protons, iVp, and therefore the ratio Nn/Np was plotted against 
i/Fp. The experimental values of ref. pO|, corrected in ref. ||121|| , give Nn/Np 



(2530 ± 1050)/(1112 ± 130). The error band in Nn/Np allowed values for F„/Fp in 
the range 0.0 — 1.65 | p.23| ], perfectly compatible with the OPE predictions. Note that 
the inclusion of the 2A^-induced channel extended the possible values for F„/Fp at 
both ends of the interval, with respect to what would have been obtained omitting 
this channel. 

The fact that the relative error of Nn/Np is the sum of the relative errors of Nn and 
Np, together with the fact that neutrons are measured with little precision, makes 
the uncertainty on Nn/Np very large and therefore leads to a poor determination 
of Tn/Tp. It was pointed out in ref. | ]123| ] that the separate measurement of the 
number of protons and neutrons per decay event, Up and n„, would provide more 
reliable predictions. As an example, fig. |TB| shows the number of protons (upper half) 
and neutrons (lower half) per decay event for different energy cut-offs (from top to 
bottom: 0, 30 and 40 MeV) in the decay of ^°^Pb. If one considers a hypothetical 
situation in which the values of Up and n„ could be measured with 10% accuracy 
and the value of F„/Fp was around 1, then this would imply having measured values 
of Up = 0.37 ± 10% and = 1.03 ± 10%, assuming an energy cut-off of 40 MeV. 
Representing these values in the figure and translating their error bands into the 
corresponding ones for F„/Fp, one would obtain, from Up, the value F„/Fp = 1.0lo;22 
and, from n„, the value F„/Fp = 1.0lo;52. This example clearly demostrates that the 
measurement of Up is, so far, the most crucial one in order to determine the ratio 
F„/Fp. An experimental set up to be used at TJNAF and at the Yerevan Electron 
Synchroton (YESS) facility has already been devised in order to measure Hp from the 



delayed fission events from the decay of heavy hypernuclei ||125|| , which hopefully will 
lead to a cleaner determination of the ratio F„/Fp. 
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Figure 16: Number of protons (n^) and neutrons (n„) per A decay event in ^°^Pb. 
Prom top to bottom the results include energy cuts of 0, 30 and 40 MeV. Dashed 
lines: IN- induced mechanism. Solid lines: lA^ + 2N mechanisms. 
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4 Conclusions and perspective 



The discussions along this review have shown some interesting features that we sum- 
marize here. 

The mesonic decay of A hypernuclei has offered some evidence of the strong re- 
pulsion at short distances of the AA^ interaction which appears naturally in current 
models of the YN interaction. The sensitivity of the mesonic decay to the pion nu- 
cleus interaction has been stressed repeatedly in the calculations. The total mesonic 
decay rate is appreciably enhanced due to the p-wave attraction of the optical po- 
tential. However, exclusive mesonic decay into particular final nuclear states may 
select the s-wave part and the partial rate can be reduced, particularly if the final 
nucleus is a closed shell one. Thus, investigations on A mesonic decay into different 
channels can offer complementary information on the pion nucleus optical potential 
and should help in deciding between different parametrizations of theoretical models, 
which are equally successful in explaining pion nuclear scattering or pionic atoms 
data. Moreover, the mesonic decay rates of heavy hypernuclei, even if small, would 
be most welcome since the pion renormalization is most important there. It should 
also be stressed that these measurements would also provide information about the 
tt'^ nuclear interaction which is not available from scattering data. 

The non-mesonic decay has been the object of intensive theoretical investigations 
and models going beyond the traditional one pion exchange mechanism are now avail- 
able. Yet, the inclusion of heavier mesons or multimeson exchange does not change 
appreciably neither the decay rate nor the ratio F^/Fp, which was the main motiva- 
tion to develop such models. The incorporation of quark degrees of freedom into the 
picture brings interesting new features, such as important A7 = 3/2 contributions to 
the decay rate and the possibility of obtaining increased values for the F„/Fp ratio. 
However, more theoretical investigations are needed to understand both the mesonic 
and the non-mesonic decays from the same basic quark Hamiltonian. Although the 
discrepancy between the theoretical models and experiment for the value of F„/Fp is 
still a puzzle, one must note that the experimental uncertainties for this ratio are very 
large. The two-nucleon induced decay was proposed as a hope to resolve this puzzle 
but, as discussed here, the consideration of this channel leads to enlarged error bars in 
the experimental analysis. Nevertheless, some efficient methods to determine F„/Fp 
with far more precision than has been possible until now have been proposed. One 
possibility is through the determination of the number of protons per A decay and 
their spectra. New experiments along these lines are now planned and most awaited. 
Another interesting feasible experiment is the measurement in coincidence of two nu- 
cleons emerging back-to-back. This should allow to disentangle the nucleons of the 
one-nucleon-induced mechanism from those of the two-nucleon-induced one, leading 
to a cleaner determination of the F^/Fp ratio. 

Measurements of the lifetime of heavy hypernuclei are now available from different 
experiments and, although agreeing at the qualitative level, more precise determina- 
tions of these rates would be also helpful to establish clearly when saturation is reached 
and, therefore, obtain information on the range of the weak AA^ — > A^A" transition. 

The A7 = 1/2 violation in the non-mesonic weak decay of hypernuclei is an 



55 



intriguing possibility that should be explored further through ratios of partial decay 
rates of light hypernuclei. 

The possibility of producing polarized hypernuclei with the {7r^,K~^) reaction 
has extended the available observables to asymmetries of the weak decay products. 
These results bring complementary information to that obtained from partial and 
total non-mesonic decay rates of hypernuclei and some theoretical models have now 
predictions for the asymmetry of the emitted particles. However, in order to constrain 
these models, it is necessary to reduce the experimental uncertainties and clarify the 
contradicting experimental results. 

In summary, the field of A hypernuclei and their weak decay modes has experienced 
an impressive progress in the last decade. The availability of new data from practically 
all intermediate energy experimental facilities is building up the grounds for a bright 
future in which one may expect to settle the present problems and open the door to 
novel an interesting physics. 
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